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ABSTRACT 


Under  a  previous  effort  "Antenna  Alignment  Study,"  F30602- 74-00276,  a 
technique  was  formulated  to  align  troposcatter  antennas  by  swinging  the  beams 
while  noting  the  symmetry  of  the  received  Doppler  spectrum.  The  measure  of 
symmetry  was  computed  from  the  cross-correlation  of  the  quadrature  components 
of  the  received  complex  envelope  signal.  A  concern  in  this  study  was  the  need 
for  synchronization  of  the  transmitter  and  receiver  clocks.  This  report 
details  the  computation  of  the  Doppler  spectrum  skew  which  is  fairly  insensitive 
to  clock  offset.  A  microprocessor  was  designed  and  programmed  to  provide  a 
real-time  skew  estimate  that  could  be  used  as  an  alignment  indicator.  Tests 
on  the  RADC  Verona-Ontario  link  employing  TRC-97  radios  showed  that  alignment 
accuracy  better  them  +  1/3®  or  approximately  +  1/10  beam-width  could  be 
accomplished  in  about  10  minutes.  It  was  found  necessary  to  set  the  average 
frequency  shift  equal  to  zero  by  adjusting  the  receiver  clock  in  order  to 
estimate  spectrum  skew  corrently  using  the  algorithm  presented  here. 
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Evaluation 


This  report  summarizes  work  aimed  at  providing  a  rapid  and  accurate  align¬ 
ment  technique  for  troposcatter  antennas.  An  Initial  program,  "Antenna  Align¬ 
ment  Study,"  performed  under  contract  F30602-74-C-0276  established  a  basis  for 
a  new  alignment  procedure  that  Is  based  on  received  spectrum  symmetry  rather 
than  received  power  level  techniques  presently  employed.  This  Initial  work 
though  was  validated  using  a  Hewlett  Packard  correlator  which  was  not  optimi¬ 
zed  for  Inclusion  Into  existing  troposcatter  radio  systems. 

The  work  under  the  present  contract  was  aimed  at  providing  for  an  alter¬ 
nate  Implementation  that  will  ultimately  provide  for  a  fleldable  piece  of 
hardware  that  will  accomplish  the  goals  of  rapid  and  accurate  alignment.  Be¬ 
sides  providing  for  rapid  alignment  It  Is  Important  to  realize  that  the  tech¬ 
nique  provides  for  accurate  alignment  on  the  great  circle  path  while  power 
level  based  techniques  can  only  provide  for  short  term  alignment  which  Is  not 
necessarily  the  great  circle  path  and  thus  requires  updating. 

Before  fielding  several  Issues  require  further  work,  these  are  validation 
of  the  alignment  technique  over  the  full  breadth  of  troposcatter  medium  varia¬ 
tions,  reworking  of  the  experimental  model  to  eliminate  erroneous  results  due 
to  aircraft  and  further  optimizing  the  Implementation. 

Based  on  results  and  the  recommendations  for  future  work  the  Improve¬ 
ment  to  troposcatter  alignment  Is  viable  and  much  closer  to  delivery  to  the 
field. 

PAUL  SIERAK 
Project  Engineer 
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Section  1 

INTRODUCTION  AND  SUMMARy 

Objectives  of  the  Program 

Until  recently,  azimuthal  alignment  of  troposcatter  antennas  was  accomplished 
by  averaging  received  power  to  find  the  proper  pointing  angle  of  each  antenna. 

The  procedure  is  arduous,  slow,  and  often  inaccurate  because  of  the  deep,  per¬ 
sistent  fading  present  in  the  channel.  In  earlier  work,  Birkemeier  and  Sill  [11 
established  a  new  technique  for  alignment.  According  to  their  method,  a 
monochromatic  signal  is  transmitted  and  the  shape  of  the  received  Doppler  spectrum 
is  used  to  determine  alignment.  The  original  method  involved  moving  the  receiving 
antennas  until  the  Doppler  spectrum  was  symmetric  about  zero  frequency.  When  this 
was  accomplished,  the  receiving  antenna  was  defocused  and  both  antennas  moved 
synchronously  clockwise  or  counter-clockwise  until  the  Doppler  spectrum  was  again 
symmetric.  The  beams  were  then  shown  to  be  pointing  along  the  great  circle  path. 

The  first  Doppler  measurements  were  obtained  with  a  mini-computer  that 
performed  an  FFT  on  the  incoming  data.  This  method  required  expensive  equipment 
and  operator  interpretation  of  the  shape  of  the  spectrum.  A  less  expensive 
method  requiring  little  operator  skill  was  sought.  The  next  approach  made  use 
of  the  fact  that  if  the  spectrum  is  not  symmetric  about  zero  frequency,  there 
will  be  a  non-zero  average  frequency.  The  average  frequency  can  be  used  as  a 
measure  of  the  symmetry  of  the  Doppler  spectrum  and  thus  of  antenna  alignment. 
There  are  several  known  algorithms  for  the  estimation  of  the  average  frequency. 
Birkemeier,  Sill,  and  Fontaine  adapted  the  correlation  method  [2]  so  that  it 
could  be  run  in  real  time  on  a  microprocessor.  However,  the  scheme  has  a 
particular  disadvantage  in  that  the  average  Doppler  frequency  will  be  in 
error  by  an  amount  equal  to  any  frequency  difference  between  the  transmitter 
and  receiver  clocks.  Thus  a  substantial  alignment  error  would  occur  with  this 
method  if  the  clocks  were  not  in  synchronism  as  shown  by  Sill  [2]. 


In  this  program  we  provide  a  microprocessor  design  together  with 
software  to  estimate  the  3rd  central  moment  or  the  skew  of  the  Doppler 


spectrum.  The  skew  provides  an  alignment  indicator  that  is  a  direct 
measure  of  the  symmetry  of  the  spectrum  about  its  mean  frequency  and 
therefore  does  not  depend  on  exact  clock  synchronization. 

In  order  to  use  the  microprocessor  with  the  TRC-97  radio  it  was  also 
necessary  to  provide  coherent  signal  sources  for  the  TRC-97  equipment. 

This  report  also  details  the  design  of  this  auxilliary  equipment. 

Technical  Approach  to  Doppler  Moment  Computations 

In  order  to  compute  the  skew,  we  will  need  to  determine  various  spectral 
moments  of  the  input  process  and  it  is  worthwhile  to  review  the  significance 
of  these  moments. 

The  first  moment,  <a,  is  the  sum  of  two  independent  frequency  terms. 

The  first  component  is  due  to  the  frequency  difference  between  the  receiver 
and  transmitter  clocks,  and  the  second  term  is  due  to  the  lack  of  symmetry 
in  the  spectrum. 

—  2 

The  second  spectral  moment  (<o-co)  is  a  measure  of  the  spread  of  the 
spectrum  about  the  mean  frequency. 

The  third  central  moment  (oi-<o)  3  is  used  to  estimate  the  symmetry  of  the 

Doppler  spectrum  about  its  mean  frequency.  If  the  moment  is  zero,  the  spectrum 

is  symmetric  about  the  mean;  if  positive,  the  spectrum  is  skewed  positively; 

if  negative,  the  spectrum  is  skewed  negatively. 

We  begin  by  defining  the  skew  to  be: 

00 

(oj-u)  3  s  (o>)  du 

Skew  =  -  ■  . . .  ■  (1) 

00 

j  S  (CO)  dw 

-CD 

where  u>  is  the  average  frequency  defined  by  [4] 


3 


£0  = 


&i  S(oi)doi 


S  (a) )  doi 


(2) 


Expanding  (1)  yields 


f  3 


] 


S  (co)  doi  -  3oj  oi  S(oi)  doi  +  35  oi  S(oi)  doi  5  S(oi)  do) 


Skew  = 


S  (01)  doi 


(3) 


We  now  note  that  the  received  signal  can  be  written  as  a  narrow  band 
process  [51 . 


S(t)  =  X(t)cos  oi  t  +  Y(t)sin  oi  t 
c  c 

joi  t  (4) 

=  Re{ [x (t)  +  jY (t) ]e  c  } 

The  quantity  X(t)  +  j <Yt)  is  known  as  the  complex  envelope  of  the  received 
signal  and  X ( t )  and  Y(t)  are  known  as  it’s  "inphase''  and  "quadrature  conponents" 
respectively.  Calling  the  conplex  envelope  Z(t);  we  now  proceed  to  determine 
the  Doppler  spectrum  of  the  correlation  functions  of  X ( t )  and  Y (t)  as  follows: 

Z(t)  =  X(t)  +  j  Y  ft)  (5) 

RZZ(T)  =  E{ [X (t)  +  jY (t) ) [X (t+T )  -  jY (t+T )} 

(6) 

=  RXX(T)  +  R^CT)  +  jRvx(X)  -  jRxy(T) 

If  we  assume  that  the  amplifiers  in  the  quadrature  channels  have  equal  gain 
and  phase  characteristics,  the  assunption  RXX(T)  =  Ryy(T)  can  be  ma<ie*  we 
further  assume  that  X(t)  and  Y(t)  are  uncorrelated  random  variables,  it  can  be 
inplied  that  Rxy(T)  =  -Ryx^  15  J.  Using  these  assunptions,  (6)  can  be  sinplified 


to  yield 
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(7) 


RZZ(T)  =  2(Rxx(T)  -  3ryx(t)) 


If  we  define  S(to)  =  F(R  (T)),  then  (8)  follows: 

LtCt 


S(to)  ej“°  dto  =  27TR  „(0)  =  4ttR  (0) 
6L  XX 


(8) 


Papoulis  [5)  has  shown  that 


V(T)  =  ~rxy(t) 


(9) 


This  result  will  be  used  to  calculate  to.  The  differentiation  relationship  of 
the  Fourier  transform  allows  us  to  write  (71: 


1 

2ff 


}toS(to)eJ  dto  =  R  (r) 


(10) 


Using  this  relationship,  it  follows  that 


jtoS(to)ej“°  dto  =  2itr‘  (0) 

oL 


(11) 


toS(«)  dto  =  —  r‘_ (0) 


We  combine  (9),  (10),  and  (11)  and  use  the  fact  that  R^  (T)  is  even  to  find  to. 
Thus 


to  =  4ttRxx(0) 


(12) 


Using  (10)  twice,  we  get  the  second  moment 


to2  Sftoje^  dto  =  -2ffR  (0) 

CUt 


(13) 


Papoulis  (SI  has  also  shown  that 


RZZ(T)  =  RZZ(T) 


(14) 


Since  RXy(T)  is  an  odd  function,  even  derivatives  at  T=0  are  zero.  This  allows 
us  to  combine  (13)  and  (14)  to  yield: 


u2  =  4ttRxx  (0) 


(15) 


Finally,  the  third  moment  is  computed  in  the  same  fashion. 


w3  S(ai)e^°  dw  =  R22(0) 


(16) 


Rz'z  <0)  *  Rzz<°> 


(17) 


W  -  4TIRxy(0) 


(18) 


We  now  combine  (8) ,  (12)  ,  (15)  ,  (18) ,  and  (3)  to  yield  the  formula  for  the 


skew: 


Skew  = 


Rxv(0) 

Rxx(0) 


3rxy(0)rxx(0) 

Rxx(0,Rxx<0) 


2(Rxt(0)l3 

tRxx(0))3 


(19) 


This  is  the  theoretical  formula  for  the  relative  skew.  Several  assunp- 
tions  about  the  process  have  been  made.  They  are  RXX<T)  "  Ryy(T),  ^Y*1*  “ 
-Ryx(T)  and  that  the  X(t)  and  Y(t)  processes  are  jointly  stationary.  For 
the  computation  algorithm  to  function  properly,  all  of  these  assumptions  must 
hold. 


Practical  Implementation  of  the  Skew  Estimator 

The  validity  of  each  of  the  assumptions  made  earlier  must  be  examined 
in  light  of  the  data.  The  most  basic  assumption  made  earlier  is  that  the 
in-phase  and  quadrature  components  of  the  signal  are  uncorrelated  random 
variables,  that  is  Rxy(0)  =  0*  It  is  clear  that  if  both  channels  possess  a 
DC  term,  this  assurption  cannot  hold.  All  amplifiers,  no  matter  how  carefully 
designed,  will  have  some  DC  offset.  This  is  especially  true  of  amplifiers  that 
must  operate  at  frequencies  close  to  zero.  The  software  must  preprocess  the 


data  to  remove  the  offset 
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The  second  major  assumption  is  that  Rxx<°)  =  K^tO).  T*1*-3  statement 
implies  that  the  statistics  of  the  channels  are  equal  and  that  the  gains  of 
the  amplifiers  in  both  channels  are  equal.  If  this  is  true,  we  should  be 
able  to  interchange  the  roles  of  X(t)  and  i  (t)  and  get  the  same  results. 

This  does  happen  with  data  that  is  not  properly  preprocessed.  Certain 

terms  are  immune  to  differences  in  the  X  and  Y  channel  gains.  In  practice, 

terms  that  contain  both  X  and  Y  are  immune  to  interchanging  roles.  Terms 
such  as  RXX(T>  or  ryy*T*  do  not  have  mixed  terms  and  will  be  vulnerable 
to  both  short  term  and  long  term  variations  in  the  gain.  These  terms  are 
particularly  sensitive  because  they  contain  derivatives  which  accent  the 
temporal  variations  in  channel  gain. 

The  long  term  variation  in  the  difference  in  the  amplifier  gains  can 
be  taken  care  of  in  software  by  correcting  the  gain  of  the  Y  channel  by 

the  ratio  of  the  rms  value  of  X ( t )  and  Y(t).  This  is  slow  age  which  cannot 

compensate  for  the  short  term  variations  in  Rxx<0)  and  R^tO).  To  correct 
for  these  short  term  variations  we  have  obtained  good  results  by  taking 
the  geometric  mean  of  Rxx<°)  and  Ryy(0)  and  usin9  this  number  to  replace 
Rxx(0)  in  the  2nd  term  in  (19). 

The  following  section  describes  the  microprocessor  implementation  of 
these  programming  objectives. 

Summary  of  Results 

Tests  of  the  Slope  algorithm  on  the  RADC  link  indicated  that  alignment 
could  be  accomplished  to  within  about  0.3°  or  about  1/10  of  a  beam  width 
with  a  period  of  about  ten  minutes.  The  Slope  Algorithm  actually  measures  <o, 
the  average  Doppler  shift  and  its  success  depends  on  exactly  synchronized 


clocks 


The  skew  algorithm  indicated  that  the  antennas  could  be  aligned  to 
1/3  degree  or  better,  also  in  a  time  period  of  about  ten  minutes.  These  measure¬ 
ments  were  based  on  single  observations  of  displayed  skew  which  occur  once  each  second. 


Visual  averaging  would  lead  to  substantial  improvement  in  expected  accuracy, 
but  such  improvement  is  not  quantified  in  this  report.  Exact  clock  synchronization 
was  not  required  in  alignment  using  skew  indication.  Instead,  w  was  set  equal 
to  zero  by  varying  the  receiver  clock  before  each  skew  measurement  was  made. 

Organization  of  Report 

The  microprocessor  hardware  and  software  design  is  summarized  in  Section  2. 
Section  3  details  the  design  of  the  frequency  synthesizers  and  detector  components 
necessary  to  interface  the  microprocessor  to  the  TRC-97  and  the  TRC-170  equipment. 
The  antenna  alignment  tests  on  the  RADC  link  are  presented  in  section  4. 

Appendix  A  contains  the  microprocessor  drawings.  Appendix  B  lists  processor 
program.  Appendix  C  presents  complete  drawings  concerned  with  the  design  and 
fabrication  of  the  synthesizer  modules  together  with  parts  lists.  Appendix  D 
contains  a  more  detailed  theoretical  description  of  the  nature  of  troposcatter 
Doppler  spectrum  and  the  estimation  of  its  parameters. 
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Section  2 

MICROPROCESSOR  HARDWARE  AND  SOFTWARE 

Introduction 

A  block  diagram  of  the  microprocessor  hardware  is  shown  in  Fig.  2.1. 

The  major  components  consist  of  the  8085  microprocessor,  4-2716  PROMS  (one 
with  test  data,  one  with  a  debugging  monitor  and  two  with  the  alignment 
programs) ,  AM9511  arithmetic  processor,  AD571  10-bit  A/D  converter,  8251  serial 
terminal  andphone  line  interfaces  for  the  program  development  and  front  panel 
switch,  light  and  BCD  display  hardware.  The  schematics  for  the  microprocessor 
system  are  in  Appendix  A.  The  system  operates  in  five  modes  (reset,  clock, 
slope,  correlation,  and  skew)  that  are  front  panel  push-button  selectable  by 
the  operator. 

The  reset  mode  resets  all  hardware  in  the  system  and  then  uses  the 
slope  algorithm  to  process  five  seconds  of  digitized  test  data  that  is  stored 
in  a  PROM.  This  provides  a  self  test  for  the  system.  The  clock  mode  is  used 
to  measure  and  display  the  offset  frequency  between  the  transmitter  and  receiver 
oscillators  while  the  slope,  correlation  and  skew  modes  are  the  three  alignment 
algorithms.  Front  panel  switches  are  scanned  at  a  ten  per  second  rate  by  the 
program  and  whenever  a  switch  is  found  pressed  the  program  is  restarted  using 
the  algorithm  that  is  indicated  by  the  switch. 

A  five  diqit  BCD  display  provides  the  operator  with  a  floating  point 
(four-digit  mantissa,  one-digit  exponent)  digital  readout  of  alignment 
information.  In  all  cases  the  displayed  quantity  represents  an  average 
that  has  been  exponentially  averaged  with  a  time  constant  of  several  seconds. 
Although  the  time  constant  is  several  seconds  the  display  is  updated  once 
a  second  with  the  new  value  of  the  running  average. 

The  arithmetic  processor  allows  arithmetic  expressions  to  be  pro¬ 
grammed  very  easily  —  much  like  a  hand  calculator.  The  10-bit  digitized 
values  for  the  in  phase  and  quadrature  signals  are  immediately  converted 


Figure 
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to  a  32-bit  floating  point  format  by  the  arithmetic  processor  and  all  sub¬ 
sequent  processing  of  these  is  done  using  the  floating  point  format. 

Two  programs  are  listed  in  Appendix  B.  One  in  FORTRAN  and  one  in  8085 
assembly,  are  included.  The  FORTRAN  program  was  used  on  a  Raytheom  704  with  tape 
recorded  data  to  develop  the  processing  algorithms.  A  deliberate  attempt  was 
made  to  maintain  a  close  correspondence  between  the  FORTRAN  and  the  8085  code 
and  in  most  cases  each  line  of  the  FORTRAN  program  is  included  as  a  comment  in 
the  8085  instructions  that  are  used  to  implement  the  FORTRAN  expression. 

General  Signal  Considerations 

The  inphase  and  quadrature  signals  are  filtered  with  analog  filters 
before  they  are  digitized  by  the  A/D  system  on  the  microprocessor  board. 

The  pass  band  of  these  filters  is  0  to  30  Hz  when  the  clock  and  skew 
algorithms  are  used  and  0.3  to  20  Hz  when  the  slope  and  correlation  algorithms 
are  used.  In  this  section  and  in  the  program  listings  the  inphase  and  quadrature 
signals  will  be  called  X  and  Y  respectively. 

The  tape  recorded  data  used  in  the  FORTRAN  program  was  sampled 
at  a  100  sample  per  second  rate.  The  programs  were  written  with  the 
intent  of  using  a  timer-drive  microprocessor  interrupt  to  sample  X  and  Y 
every  0.01  seconds ,  accumulate  data  from  100  sarrples,  and  update  the  display 
once  a  second  with  the  new  data  included  in  the  running  average.  The  timer 
interrupt  was  not  implemented  because  the  microprocessor  program  when  operating 
at  full  speed  result  in  approximately  a  100  sample  per  second  sampling  rate. 

The  displayed  variables  for  clock  or  antenna  alignment  algorithms  are 
averaged  quantities  and  hence  statistical  in  nature.  Two  conditions  exist 
that  can  introduce  large  errors  in  these  averages.  One  situation  is  when  X 


and  V  simultaneously  become  very  small  and  the  other  is  reflections  from 
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aircraft.  The  first  situation  is  easily  programmed:  the  incoming  data  are 
ignored  whenever  the  magnitude  of  both  X  and  Y  are  less  than  2%  of  the  maximum 
value. 


The  aircraft  contamination  is  more  difficult  to  detect.  Aircraft  introduce 
higher  than  normal  frequencies  into  X  and  Y  and  these  are  particularly  trouble¬ 
some  in  the  skew  algorithm  that  requires  the  second  derivatives  of  X  and  Y. 

In  an  attempt  to  eliminate  data  with  aircraft  contamination  the  frequency  X 
and  Y  is  estimated  from  each  successive  pair  of  sanfiles  and  the  data  ignored 
whenever  the  frequency  exceeds  8  Hz.  The  first  and  second  derivatives  are 
calculated  using 


Axn  = 


Vi  -  Vi 


A  XN  =  Vl  "  2XN  +  Vl 


Thus  whenever  data  are  ignored,  three  valid  samples  are  required  before  the 
derivative  calculations  can  continue. 

Measurements  made  on  the  inphase  and  quadrature  detector  outputs  indicated 
that  both  contained  non-zero  dc  components  and  that  the  gain  of  the  two  detectors 
varied  slightly.  Both  effects  are  removed  in  the  processing  by  computing  long 
time  dc  and  mean  square  values  of  the  detector  outputs  and  equalizing  the 
sampled  values  with  them: 


As  mentioned  earlier  an  exponential  averaging  technique  is  used  to 

~2  • 

form  a  running  average  of  x  t  skew,  etc.  If  A^_^  is  the  current 

value  of  the  running  average  and  I„  is  the  new  input  to  be  added  to  the 

N 

average  the  new  average  is  computed  from 


=  A. 


XN  - 


N-l 


N-l 
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8085  Processing  -  Overall  View 

The  four  FORTRAN  subroutines  described  above  perform  identical  processing 
in  the  8085  microprocessor  system.  A  brief  description  of  additional  routines 
is  given  in  this  section.  The  X,  Y  sanple  pair  is  obtained  by  the  subroutine 
ATOD.  The  routine  will  return  either  test  data  that  is  stored  in  PROM  or  will 
sample  live  data  with  the  A/D  system.  Test  data  is  returned  if  the  AC  power 
system  has  just  been  turned  on  or  if  the  operator  has  pushed  the  RESET  button; 
otherwise,  sampled  live  data  is  returned. 

The  front  panel  switches  are  scanned  with  the  routine  RDFRNT.  If  a  switch 
is  found  depressed  the  system  is  reset  and  a  logical  flag  corresponding  to  the 
switch  is  set.  Other  subroutines  use  the  flags  to  determine  the  processing  that 
is  required. 

A  scheduler  routine,  SCHED,  is  used  to  direct  the  real-time  processing. 
Every  call  to  SCHED  results  in  a  call  to  ATOD  for  a  new  data  sample  and  a  call 
to  PROCES  to  process  the  data  sanple.  Every  10th  call  to  SCHED  results  in  a 
call  to  RDFRNT  to  scan  the  front  panel  switches  and  every  100th  call  produces 
a  call  to  UPDATE  and  DISPLY.  As  mentioned  earlier  the  system  was  designed 
for  a  programmable  timer  to  produce  a  call  to  SCHED  every  .01  sec  but  this 
is  not  implemented.  Instead  each  return  from  SCHED  is  immediately  followed  by 
another  call  to  this  routine. 

Program  Implementation  of  Clock  and  Alignment  Algorithms 

Under  ideal  conditions  of  noise  free,  zero  doppler  shift,  transmission 
with  the  transmitter  and  receiver  oscillators  operating  at  a  frequency  dif¬ 
ference  of  oj.,  the  X  and  Y  signals  become: 
d 


x  (t )  =  cos  o>dt  and  y(t)  =  sin  M^t 
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where  N  is  the  number  of  samples  that  have  been  taken.  This  algorithm  acts 
like  an  RC  averaging  circuit  with  a  time  constant  of  NT  seconds  where  T  is 
the  time  between  new  inputs  1^.  Thus  new  information  is  weighted  heavily 
and  old  information  is  forgotten.  In  addition  N  is  varied  during  the  measure¬ 
ment  starting  at  a  value  of  one  and  increasing  to  a  maximum  of  16  or  60 
depending  on  the  variable  that  is  averaged, 

FORTRAN  Processing  -  Overall  View 

The  FORTRAN  program  has  four  subroutines  INITAL,  PROCES,  UPDATE,  and 
DISPLY  for  processing  the  data  and  each  of  these  has  a  counterpart  in  the 
8085  program.  A  program  listing  of  the  FORTRAN  program  is  in  Appendix  A. 

The  data  used  in  the  FORTRAN  is  read  in  from  magnetic  tape  in  blocks  of 
1024  digitized  (12-bit),  pairs  of  samples  X^,  by  the  main  program.  The 
main  program  then  passes  one  pair  of  samples  at  a  time  to  the  subroutine 
PROCES  just  as  the  A/D  subroutine  in  the  8085  system  would  pass  one  pair 
to  the  8085  PROCES  routine.  Each  subroutine  contains  the  processing  necessary 
for  all  four  algorithms:  clock,  slope,  correlation,  and  skew  but  uses  logical 
flags  to  select  only  the  processing  that  is  necessary  for  the  algorithms  in 
use  when  the  subroutine  is  called. 

The  subroutine  INITAL  is  used  to  initialize  all  arrays,  exponential 
averaging  parameters,  etc.  each  time  a  run  is  started  or  whenever  the  operator 
decides  to  reset  or  clear  the  program  of  accumulated  data. 

The  PROCES  routine  is  called  once  for  each  X,  Y  pair  of  samples.  Data 
are  tested  for  validity,  derivatives  calculated,  one-second  sums  of  dc  and 
mean  square  values  of  X  and  Y  are  tabulated  for  the  100  samples  that  are 
passed  to  PROCES.  After  100  samples  (equivalent  to  1  second  of  data)  are 
processed,  the  accumulated  runs  are  used  by  UPDATE  to  add  to  the  running  averages 
of  these  variables  and  the  parameters  in  PROCES  are  reset  to  allow  for  the 
accumulation  of  the  next  100  samples.  After  the  running  averages  are  updated 


the  routine  DISPLY  is  called  to  display  the  new  values 


The  frequency  can  be  obtained  by  defining  9  =  tan  ^  (y/x)  and  calculating 


dO  dy  dx 

3F  =  x  5F  '  y  dT  =  “ 


Jd  • 


In  terras  of  signals  that  are  sampled  every  At  sec  the  difference  frequency  is 
estimated  as 

A0.  x.Ay.  -  y.Ax. 

.  _  i  11  ii 


d  2irAt 


2irAt 


A0  . 


The  running  average  (16  sec  time  constant)  of  the  quantity  — is  continuously 


displayed  when  the  clock  algorithm  is  used. 


(Ax.  )Y. 


The  slope  algorithm  displays  a  running  average  of  the  quantity  — = -  . 

x2At 

The  X^  and  Y^  used  in  this  expression  and  in  all  calculations  are  the  sample  values 

— 2 

with  the  dc  removed  and  normalized  to  the  same  mean  square  value.  The  x  term  is 
a  running  average  of  the  mean  square  value  of  x  that  is  used  to  normalize  the 
slope  statistic  so  that  it  is  independent  of  signal  strength. 

The  cross  correlation  between  continuous  function  of  time  x(t)  and  y(t) 


is  defined  as  R  (T)  =  x(t)y(t-T).  For  sanpled  signals  the  cross  correlation 
xy 

becomes 


R^trr^iT}  =  x  (kAt)y  (kAt-mAT)  . 

The  cross  correlation  is  computed  for  -1.5  <  T  <  1.5  seconds  by  using 

At  =  At  =  0.03  sec  and  a  table  to  store  the  100  values  for  R  (mAx)  that 

xy 

correspond  to  -50  <  m  <  50.  Since  .03  is  three  times  the  sampling  interval, 
every  third  sample  of  raw  data  is  used  in  calculating  the  cross  correlation. 
These  data  samples  are  placed  in  a  "delay  line",  first  in  first  out  type 
table,  50  samples  long  for  X  and  100  samples  long  for  Y.  As  each  X  sample 
reaches  the  bottom  of  the  table  it  is  correlated  with  the  50  samples  of  Y  for¬ 


ward  and  backward  in  time  and  then  discarded.  Since  the  alignment  criteria 
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with  the  correlation  algorithm  is  to  minimize  the  absolute  slope  about  T=0, 

the  slope  of  a  line  from  the  origin  to  the  maximum  value  of  R  (mAx)  is  the 

xy 

quantity  that  is  displayed  to  the  operator. 

The  equation  for  skew  was  developed  earlier;  and  although  it  is  a 
straightforward  arithmetic  calculation,  it  is  a  sensitive  measure  with  a 
larger  variance  than  the  slope  or  correlation  measures.  This  is  due  to  the 
subtraction  of  terms  in  the  expression  that  frequently  are  very  large  and  nearly 
equal. 

Summary  and  Recommendations 

The  hardware  and  software  have  proven  to  be  very  trouble  free  and  reliable. 
The  higher  doppler  frequencies  that  were  present  in  data  on  the  RADC  link,  but 
not  on  the  UW  link,  frequently  caused  the  system  to  reject  the  data  as  aircraft 
contamination.  Other  means  of  detecting  aircraft  reflections  should  be  tried 
and  the  most  promising  of  these  is  to  look  for  rapid  fluctuations  in  envelope 
that  exceed  a  peak  to  peak  threshold. 


INTERFACING  WITH  TRC-97  AND  TRC-170 


Frequency  Synthesizer  Design 

The  Air  Force  troposcatter  TRC-97  and  TRC-170  equipment  requires  a 
frequency  synthesizer  that  will  convert  a  precision  frequency  reference  into 
coherent  70  MHz  and  50  MHz  outputs.  The  receiver  also  requires  a  quadrature 
detector  for  its  70  MHz  intermediate  frequency  signal.  A  block  diagram  of 
the  system  is  shown  in  Appendix  C,  Figures  Cl  and  C2.  As  the  figures  show 
t  'e  radios  are  to  use  their  own  internal  50  MHz  and  70  MHz  sources  when 
operated  in  the  non-alighment  mode  and  the  synthesizer's  coherent  oscillator 
signals  in  the  alignment  mode.  Coaxial  relays  are  included  in  the  aligner  for 
switchiig  these  signals.  The  50  MHz  and  70  MHz  synthesizer  outputs  must  have 
a  frequency  stability  of  at  least  1  part  in  10  ,  and  both  the  50  MHz  and 

the  70  MHz  outputs  must  provide  18  dbm  to  a  50  ohm  resistive  load.  Harmonic 
content  within  10  MHz  of  the  fundamental  was  to  be  down  at  least  20  db. 

These  specifications  were  met  with  an  "indirect"  frequency  synthesizer, 
employing  phase  locked  loops  (PLL)  for  frequency  multiplication.  The  basic 
designs  for  these  synthesizers  are  shown  in  Figure  C3  and  C4. 

The  Precision  Quartz  Clock 

The  precision  5  MHz  crystal  oscillator  which  provides  the  source  for  the 
phase  lock  loops  is  a  double-oven  Vectron  model  CO-246V-1  oscillator  with 
frequency  stability  of  1  part  in  10^.  The  oscillator  frequency  cam  be 
externally  controlled  by  a  DC  voltage.  The  connection  diagram  for  the 
frequency  standard  is  shown  in  Figure  C5. 

The  25  MHz  and  35  MHz  Phase  Lock  Loops 

The  5  MHz  reference  signal  is  power  divided  between  the  50  MHz  and 
70  MHz  phase  locked  loops.  The  basic  phase  locked  loop  design  is  shown  in 
Figure  C6.  The  voltage  controlled  oscillators,  VCO's,  oscillate  at  half 
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the  desired  output  frequency  (35  MHz  for  the  70  MHz  unit  and  25  MHz  for  the 
50  MHz  unit)  because  integrated  circuit  PLL's  are  limited  to  a  maximum 
frequency  of  about  35  MHz .  A  balanced  mixer  is  used  to  double  the  VCO  frequencies 
to  the  desired  70  MHz  or  50  MHz  values. 

The  XR215  PLL  integrated  circuit  chosen  consists  of  a  phase  detector  and 
VCO.  The  lowpass  filter  corner  frequency  for  the  control  signal  and  the  VCO 
frequency  range  are  determined  by  an  external  capacitance  connected  to  the 
XR215  chip.  The  frequency  divider  is  the  93S05  variable  modulo  counter.  For 
the  70  MHz  circuit  it  is  wired  to  frequency  divide  by  7;  for  the  50  MHz  circuit 
it  is  wired  to  divide  by  5. 

The  70  MHz  phase  locked  loop  schematic  is  shown  in  figure  C7.  Due  to 
the  variation  in  characteristics  from  one  VCO  to  another  the  VCO  capacitor 
C107  must  be  selected  to  give  the  correct  VCO  frequency  range.  C107  should 
be  no  larger  than  12  pF.  C105  and  C106  determine  the  comer  frequency  of  the 
tuning  voltage  low  pass  filter.  The  filter  transfer  function  has  a  single  pole 
at  1.33  x  iou  or  212  KHz.  R100 ,  the  VCO  free  run  frequency  adjusting 
resistor,  should  be  adjusted  to  allow  the  PLL  to  achieve  lock.  Q100  is  part  of  a 
of  a  common  emitter  amplifier  that  raises  the  VCO  signal  amplitude  to  the  level 
required  to  drive  the  Schmidt  Trigger,  U101.  The  Schmidt  Trigger  converts 
the  analog  35  MHz  signal  to  a  TTL  compatible  signal  to  prevent  erratic  operation 
of  the  frequency  divider  which  follows  it.  The  frequency  divider,  U102 ,  takes 
the  35  MHz  signal  from  the  Schmidt  trigger  and  divides  its  frequency  by  7. 

Thus  the  output  of  U102  has  a  frequency  of  5  MHz.  This  signal  is  coupled  to 
the  phase  detector  in  U100  where  it  is  phase  compared  with  the  5  MHz  frequency 
standard.  The  output  of  the  Schmidt  trigger  is  also  coupled  to  [101  which  is 
part  of  an  emitter  follower  power  amplifier.  This  power  amplifier  drives  a 
double  balanced  mixer,  U103,  with  the  35  MHz  signal.  The  mixer  is  used  as  a 
frequency  doubler  to  generate  the  desired  70  MHz  signal.  This  signal  is  passed 
to  the  70  MHz  amplifier  for  amplification  and  band  pass  filtering. 


The  50  MHz  phase  locked  loop  schematic  diagram  is  shown  in  Figure  t'8 , 
and  its  parts  list  follows.  The  50  MHz  PLL  is  similar  to  the  70  MHz  PLL. 

The  only  differences  are  a  larger  coupling  capacitance,  C214,  a  different 
wiring  of  the  frequency  divider  to  divide  by  5  instead  of  7  and  a  resistor, 
R218,  in  series  with  the  50  MHz  output.  This  resistor  attenuates  the  50  MHz 
output  to  prevent  it  from  saturating  the  50  MHz  amplifier  and  generating 
excessive  harmonic  levels. 

The  70  MHz  Amplifier 

The  oriqinal  design  with  minor  modifications  was  used  for  the  70  MHz 

amplifier.  The  schematic  diagram  is  shown  in  Figure  C9.  The  output  of  the 

amplifier  goes  to  the  detector  in  the  receiver  and  to  jack  J1  on  the  trans¬ 
mitter  A-l 7  assembly  by  way  of  the  70  MHz  coax  relay. 

The  original  design  with  minor  modification  was  also  used  for  the  50  MHz 

amplifier.  The  schematic  diagram  is  shown  in  Figure  CIO.  R400  was  added  to  the 

amplifier  input  to  eliminate  self  oscillation.  L400  and  C400  were  also  added 
to  the  amplifier  input  as  a  series  tuned  circuit  to  reduce  the  75  MHz  spurious 
output  from  the  50  MHz  phase  locked  loop. 

The  Quadrature  Detector  and  Low-Pass  Filter 

The  receiver  also  requires  a  synchronous  quadrature  detector.  This 
detector  must  demodulate  the  70  MHz  intermediate  frequency  signal  into  its 
quadrature  components.  These  signals  must  then  be  low  pass  filtered  with  a 
comer  frequency  of  20  Hz.  The  schematic  diagram  of  the  detector  is  shown 
in  Figure  Cll.  The  low  pass  filter  contained  in  the  detector  was  designed 
to  have  a  voltage  gain  of  10.  This  should  result  in  2  volt  peak-to-peak  output 
levels  from  the  X  and  Y  channels  under  nominal  channel  conditions.  This  is 
the  level  required  to  drive  the  following  stage  which  is  a  DC  blocking  filter. 
The  low  pass  filter  is  a  Sallen  and  Key  configuration  designed  to  have  the 
following  transfer  function: 


H(s) 
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This  results  in  a  two  pole  low  pass  filter  with  a  comer  frequency  of  20  Hz 
and  a  Q  of  0.707.  This  Q  gives  the  steepest  possible  skirt  without  any 
passband  ripple. 


Low-Pass,  Band-Pass  Filter 

The  output  of  the  detector  module  is  next  fed  to  the  Low-Pass,  Band-Pass 
Filter.  One  section  of  this  filter  additionally  attenuates  frequencies  above 
20  Hz  by  6  db/octave  and  outputs  this  low-pass  signal  to  the  microprocessor 
to  be  used  in  the  "slope"  program  to  compute  the  average  Doppler  frequency 
or  the  slope  of  the  cross-correlation  function  at  zero  lag.  The  band-pass 
filter  section  adds  a  high  pass  filter  in  cascade  to  the  low-pass  section  to 
remove  frequency  components  below  0.3  Hz.  The  band-pass  filtered  outputs  are 
fed  separately  to  the  microprocessor  to  be  used  to  compute  the  average 
frequency  whenever  a  strong  specular  signal  at  zero  Doppler  shift  exists  in 
the  signal.  In  tests  at  RADC  this  output  was  not  used.  The  schematic 
diagrams  and  parts  layout  are  found  in  Fig.  ell  and  Fig.  C23,  respectively. 


The  Power  Supply 

The  frequency  synthesizer  power  requirements  are  5  volts  for  the  phase 
locked  loops,  12  volts  for  the  phase  locked  loops,  amplifiers  and  detector, 

24  volts  for  the  frequency  standard  and  -12  volts  for  the  detector.  The  power 
supply  schematic  diagram  is  shown  in  Figure  C12.  The  power  supply  module 
contains  the  bridge  rectifiers,  filter  capacitors  and  24  volt  regulator.  The 
outputs  of  the  module  are  24  volts  regulated  for  the  frequency  standard  and 
relay  coils,  40  volts  unregulated  for  the  12  volt  and  5  volt  regulators  and 
-20  volts  unregulated  for  the  -12  volt  regulator.  The  12  volt,  5  volt,  and 
-12  volt  regulators  are  preceeded  by  voltage  dropping  resistors  to  reduce 
the  voltage  regulator  dissipation.  The  12  volt,  5  volt,  and  -12  volt 
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regulators  are  mounted  on  the  frequency  synthesizer  base  plate  for  heat 
sinking.  Because  the  transmitter  does  not  have  a  detector  circuit,  components 
C60 1 ,  R602 ,  CR601  and  U603  are  not  used  in  the  transmitter  power  supply. 

Each  phase  locked  loop,  amplifier  and  detector  circuit  was  built  on  its 
own  printed  circuit  board  and  mounted  in  its  own  box.  All  power  supply 
connections  to  the  boxes  are  made  with  feed  through  capacitors.  All  other 
connections  to  or  between  boxes  are  made  with  TNC  connectors  or  feedthrough 
terminals.  This  connection  scheme  is  used  to  eliminate  electro-magnetic 
interference  between  the  frequency  synthesizer  and  other  system  components. 

Test  Results 

When  the  synthesizers  were  tested,  the  phase  locked  loop  selected 
capacitors  and  VCO  tuning  potentiometers  were  adjusted  while  observing  the 
50  MHz  and  70  MHz  signal  outputs  with  a  spectrum  analyzer.  These  components 
were  adjusted  to  values  that  allowed  the  phase  locked  loops  to  achieve  lock 
at  50  MHz  and  70  MHz  after  the  internal  temperatures  stabilized  with  an 
ambient  temperature  of  25°  C.  The  50  MHz  and  70  MHz  output  spectral  purities 
and  signal  levels  are  recorded  in  photographs  C13  and  C14.  The  units  were 
each  driving  a  50  ohm  resistive  load  for  these  results.  It  can  be  seen  from 
the  photographs  that  the  signal  output  level  requirements  are  met  with  these 
synthesizers. 

Additional  Suggestions 

The  spectral  purity  of  the  frequency  synthesizers  could  be  improved 
at  an  added  cost.  The  noise  sidebands  on  the  50  MHz  and  70  MHz  signals  could 
be  improved  by  replacing  the  RC  VCO’s  in  the  phase  locked  loops  with  voltage 
controlled  crystal  oscillators,  VCXO's.  VCXO  frequencies  would  be  more  stable 
with  age  and  temperature  variations.  The  use  of  VCXO’s  should  also  reduce 
the  VCO  frequency  variations  from  IC  to  IC.  This  would  eliminate  the  need 
for  selected  capacitors  in  the  PLL  circuits.  The  XR215  can  be  used  as  a  VCXO. 
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The  spectral  components  that  are  spaced  every  5  MHz  about  the  50  MHz 
and  70  MHz  desired  signals  could  be  greatly  attenuated  by  installing  crystal 
filters  between  the  phase  locked  loops  and  amplifiers.  An  additional  stage 
would  need  to  be  added  to  compensate  the  insertion  loss. 

Printed  Circuit  Boards 

Art  work  for  PC  boards  used  in  the  synthesizers  are  shown  in  Appendix  C 
Figures  C 1S-C18. 

Parts  Layout 

Layout  of  parts  for  the  various  synthesizer  modules  are  shown  in 
Figures  C19-C22.  The  50  MHz  PLL  is  identical  to  the  70  MHz  parts  layout 
shown,  except  its  parts  are  numbered  in  the  two-hundreds. 

Parts  Lists  for  Synthesizer  and  Filter  Modules 
See  pages  C.28  through  C.36,  Appendix  C. 

TRC-170  Interfacing  Requirement 

The  TRC-170  equipment  is  fully  coherent  using  automatically  controlled 
clocks  to  drive  its  frequency  synthesizers.  As  such  there  is  no  problem  in 
interfacing  this  radio  with  the  aligner's  microprocessor.  The  aligner's 
detector  package,  however,  must  be  driven  by  the  TRC  170's  detector,  which 
in  turn  must  be  driven  by  a  70  MHz  source.  This  source  is  provided  by  the 
aligner's  synthesizer  using  its  own  internal  5  MHz  standard  or  by  substituting 
the  TRC-170 ' s  5  MHz  standard.  The  latter  would  be  preferred,  since  it 
eliminates  the  need  to  set  the  clock  according  to  skew  estimation  procedure 
outlined  in  Section  4.8. 
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Sect  ion  4 
TEST  RESULTS 

UW  Link  Tests 

The  antenna  aliqner  tests  were  performed  in  accordance  with  the  Proposed 
Test  Plan/Procedures  (clin  item  0004).  Test  efforts  on  the  940  MHz  University 
of  Wisconsin  link  held  on  12/14/79  and  12/19/79  were  largely  unsuccessful  as  a 
result  of  more-or-less  continuous  aircraft  intrusions  on  both  days.  Moreover, 
it  proved  impossible  to  skew  the  Dopple  spectrum  by  azimuthally  displacing  the 
receiving  antenna.  This  is  possibly  explained  by  the  occurrence  of  equal  and 
opposite  cross-path  wind  fields  in  scattering  layers  of  equal  cross-section, 
although  we  have  no  meteorological  proof  of  such  a  situation  existing. 

UW  Test  of  Skew  Algorithm 

We  observed  the  skew  algorithm  to  correctly  reflect  the  chanqinq  slope 
of  the  Doppler  spectrum  when  an  aircraft  crossed  the  common  volume,  but 
unfortunately,  the  skew  measure  also  was  influenced  by  the  clock  setting 
contrary  to  theoretical  expectations.  We  speculated  that  this  was  the  result 
of  a  program  deficiency  that  causes  the  1st  and  2nd  derivatives  of  the 
quadrature  components  of  the  received  data  to  be  underestimated  when  the 
fading  is  too  slow.  Fading  on  these  days  was  extremely  slow  due  to  a 
strong  specular  component  in  the  received  signal.  We  expect  faster  fading 
at  4  GHz  to  mitigate  this  problem. 

UW  Test  of  Slope  Algorithm 

Aircraft  traversals  provided  the  only  way  to  test  the  slope  algorithm. 
The  slope  changed  sign  as  the  aircraft  crossed  the  qreat-circle  in  every 
case  qualitatively  indicating  proper  operation. 
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UW  Test  of  Clock  Setting  Algorithm 

This  algorithm  was  tested  by  varying  the  5  MHz  clock  frequency  in  the 
receiver.  The  average  frequency  displayed  agreed  with  the  amount  of  observed 
clock  offset. 

UW  Test  of  Correlation  Algorithm 

The  correlation  algorithm  is  protected  from  false  data  during  on 
aircraft  intrusion  by  a  scheme  that  rejects  an  entire  30  second  data 
stream  if  the  phase  of  the  received  signal  changes  by  more  than  eight 
degrees  from  one  sample  to  the  next  (0.1  sec).  Since  aircraft  were 
particularly  numerous,  this  algorithm  succeeded  only  on  brief  occasions 
and  could  not  be  evaluated  realistically.  Later  it  turned  out  that  fast 
fading  on  the  RACD  link  was  able  to  activate  the  protective  device  and 
therefore  defeat  the  estimation  algorithm. 

Conclusions  from  UW  Tests 

While  these  tests  were  far  from  satisfactory  they  indicated  that 
we  could  expect  proper  operation  under  the  conditions  to  be  expected 
on  the  shorter,  higher  frequency  RADC  link.  We  therefore  begin  tests 
in  Verona,  Now  York  on  1/17/79. 

RADC  Link  Tests 

The  RADC  TRC-97  link  from  Ontario,  New  York  to  Verona,  New  York  was 
operated  at  4880  MHz  using  Cesiuum-beam  clocks  and  HP  synthesizers  to 
establish  all  signal  sources  in  the  transmitter  and  receivers.  Detection 
was  accomplished  using  the  coherent  quadrature  detector  packaqe  of  the 
Antenna  Aligner  together  with  the  70  MHz  reference  signals  established  by 
the  5  MHz  Vectron  oscillator  also  contained  in  the  aligner.  Observation 
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of  the  X-Y  display  of  the  quadrature  components  on  an  oscilloscope  revealed 
random  phase  jitter  of  approximately  +  45°  on  the  received  signal.  This 
was  pure  phase  jitter  that  apparently  was  generated  in  the  TRC-97  transmitter. 
(No  phase  jitter  was  observed  when  the  receiver  was  fed  the  4880  output  from 
a  coherent  source) .  This  phase  noise  caused  us  no  problem  because  most  of 
its  frequency  components  were  above  20  Hz  ar.d  were  therefore  eliminated  by 
the  low-pass  filters  in  the  aligner.  Nevertheless,  such  jitter  constitutes 
a  potential  source  for  error  if  the  filter  bandwidths  are  expanded  as 
explained  later.  Other  troubles  in  getting  the  link  up  were  due  to  cold 
environmental  temperatures  (<  40°  F)  at  the  Verona  site.  In  particular 
the  aligner's  phase-lock  loops  were  reluctant  to  lock  and  the  microprocessor 
clock  chip  refused  to  operate  until  raised  to  room  temperature. 

Test  of  Slope  Algorithm,  Step  1 

The  s  lope  algorithm  measures  the  presence  or  absence  of  a  cross¬ 
correlation  function  by  measuring  its  slope  at  the  origin.  Zero  slope  implies 
that  the  Doppler  spectrum  is  symmetrical  about  zero  frequency  shift.  In  order 
to  provide  unambiguous  results,  the  shift  due  to  clock  error  must  be  zero.  In 
these  tests  the  signal  contained  no  zero-Doppler-shi fted  component  so  the 
receiver  clock  could  not  be  set  to  zero  beat.  We  therefore  set  the  aligner's 
clock  to  zero  beat  with  the  Cesium-beam  standard  supplying  the  receiver.  The 
slope  display  was  sampled  and  recorded  each  second  for  60  seconds.  Three 
such  runs  were  made  —  the  first  two  about  one  hour  apart,  and  the  third  on 
the  next  day.  The  data  summary  is  presented  in  Table  1  below.  In  each 
case  the  transmitting  antenna  was  set  on  the  great-circle  and  both  antennas 
were  focused.  The  slope,  S,  in  volts2/sec  is  given  for  the  receiving  antenna 
displaced  +  1°  from  the  great-circle. 
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Table  1.  Slope  F.stimation  for  Stop  1  of  Alignment  Process 

Roc.  Ant.  =  1°  N  of  great-circle 

Rec.  Amt.  =  1°  S  of 

Great  Circle 

Run  1 

3:00  EST 
12/21/79 

S=-l .68 

V-27 

S/Os= 6.2 

S= 1 . 09 

V-14 

S/Os=7.8 

Run  2 

4:30  EST 
2/21/79 

S=-1.08 

Os=.33 

S/0  =  3. 20 

S=2 . 56 

Os=.22 

S/0  S=H.  6 

Run  3 

8:30  AM 
12/22/79 

l/3| 

II 

i— • 

O 

O 

CD 

Os=.33 

S/0s=  3.05 

5= 1 . 46 

Os=.39 

S/Os=3.75 

From  the  table  above  we  estimate  the  al'iqnment  accuracy  for  the  first  step 
in  the  procedure  by  noting  that  the  ratio  of  the  mean  slope  S  to  the  standard 
deviation  of  the  data,  a,.,  has  a  minimum  value  of  3.05  in  Run  3.  This  means 
that  1  degree  of  receiving  antenna  offset  resulted  in  a  mean  slope  reading 
about  three  times  larger  than  the  rms  fluctuations  in  the  displayed  slope 
data.  The  expected  accuracy  is  thus  about  1/3  of  a  degree.  That  is,  treating 
the  1  second  samples  as  independent  data,  a  single  sample  of  the  slope  measure 
has  probability  of  .64  of  lying  within  +  one  standard  deviation  of  the  true 
mean.  Note  that  in  Run  1  on  the  previous  day  about  twice  this  accuracy 
would  have  been  obtained.  Note,  however,  that  the  increase  appears  to  be  due 
to  increased  S  when  the  receive  antenna  was  1°  North  and  due  to  decreased 
rms  fluctuations  when  the  antenna  was  1°  South.  Presumably  spatial  and  temporal 
nonstationarities  in  the  scattering  structure  were  responsible  for  these 
variations.  The  minimum  obtained  accuracy  of  1/3°  is  less  than  the  back-lash 
in  the  TRC-97  antenna  gearing  and  is  therefore  probably  sufficient  in  this 
case.  The  actual  obtainable  accuracy  is  much  better,  since  the  operator 
of  the  aligner  mentally  averaqes  as  he  observes  the  displayed  numbers. 

The  rms  value  of  the  fluctuations  droo  with  1/»N  where  N  is  the  number  of 
independent  observations  averaged.  While  the  one-second  data  values  displayed 
are  not  entirely  independent,  a  substantial  improvement  is  certain  to  be 


obtained  in  actual  practice  by  a  experienced  operator. 
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Test  of  Slope  Algorithm,  Step  2 

Step  one  of  the  alignment  procedure  accomplishes  straddling  antenna 
beams.  In  step  two  of  the  procedure  the  receiving  antenna  is  defocused  and 
the  beams  are  synchronously  swung  until  the  cross-correlation  slope  is  again 
averaging  zero.  The  beams  are  then  in  great-circle  alignment.  Several 
runs  were  made  using  about  60  samples  of  slope  data  with  the  beams  displaced 
by  various  equal  and  opposite  incremental  angles  starting  at  a  maximum  of 
+  1  1/2°.  The  results  are  tabulated  in  Table  2  below.  All  five  runs  were 
made  within  one  hour  from  12:10  FST  to  13:10  EST  on  12/21/79. 


■ 

Table  2 

.  Slope  Estimation 

for  Step  2  of  Alignment  Procedure 

B 

Antenna  Setting  from 
Calibration  on  Mast. 

S 

°S 

s/os 

Rema  rks 

i 

R=  1  1/2’  N 
Focused 

T  =  1  1/2°  S 
Focused 

<3* 

O 

1 

.  308 

— 

Both  antennas  focused 

2 

T -  1 1/2°  S 
Focused 

+  1.19 

.22 

5.4 

Receiving  antenna 
de  focused 

3 

R  =  1"  N 

De  focused 

T  =  1°  S 
Focused 

+  .98 

.20 

4.9 

4 

R  =  1/2°  N 

De  focused 

T  =  1/2°  S 
Focused 

+  .11 

.24 

— 

Both  antennas  very  close 
to  great  circle 

5 

R  =  0° 

De  focused 

T  =  0° 
Focused 

-.783 

.32 

— 

Note  from  Run  4  in  Table  2  that  the  data  implies  the  correct  antenna  settings 
would  be  about  1/3°  N  for  receiving  antenna  and  1/3°  S  for  the  transmitting 
antenna.  This  is  just  the  amount  of  the  apparent  back-lash  in  the  antenna 
pointing  gear  train  and  accounts  for  the  displacement  due  to  the  wind  loading 
on  each  antenna.  (A  north  to  south  wind  prevailed  during  the  experiment.) 


rt  was  noted  in  preliminary  checks  of  the  skew  estimation  algorithm  that 


the  measure  was  sensitive  to  the  receiver  clock  setting  which  is  contrary  to 
theoretical  expectations.  It  was  possible  to  align  the  antennas  with  an 
improper  initial  clock  setting  using  the  skew  algorithm,  however.  To  accomplish 
this  we  simply  set  the  clock  to  force  the  average  "Doppler"  frequency  u)  equal 
to  zero  as  indicated  by  the  clock  mode  display  of  the  aligner.  With  <o  equal 
to  zero  the  skew  mode  of  the  aliqner  then  displays  the  true  skew  measure  of 
the  received  spectru.  To  estimate  the  alignment  accuracy  we  measured  skew 
with  the  receive  antenna  set  off  the  great  circle.  Data  for  two  antenna 
positions  are  shown  in  Table  3.  The  transmit  antenna  was  set  to  the  zero 
degree  mark  or  a  true  position  of  about  1/3°  N .  (The  wind  now  was  south  or  north 
across  the  link).  The  receive  antenna  was  set  'i  1°  S  or  a  true  position  of 
about  2/3°  S  for  Run  1.  For  Run  2  the  receive  antenna  was  set  at  1“  N  or  a 
true  position  of  about  1  1/3°  N.  The  skew  measures  correctly  reflect  these 
antenna  positions  as  shown  in  Table  3  below. 


Table  3.  Skew  Estimation  for  Step  1  of  Alignment  Procedure 

Run 

Antenna  Setting 

Both  Antenna.;  Focused 

SK 

°SK 

sk/0sk 

1 

R  =  2/3°  S 
true 

T  =  1/3°  N 
true 

-176 

H  . . 

3.95 

4.4 

2 

R  -  1  1/3°  N 
true 

T  =  1/3°  N 

true 

+428 

114 

3.75 

In  order  to  estimate  the  alignment  accuracy  from  this  data  note  from  Run  1  that 
a  displacement  of  the  receiving  antennas  of  only  1/3°  from  its  correct  straddle 
point  produced  a  significant  average  skew  value  SK  that  was  4.4  times  larger 
than  its  standard  deviation.  This  implies  an  alignment  accuracy  of  about 
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1/12  degree.  In  Run  2  the  receive  antenna  was  displaced  twice  as  far  as  in 
Run  1  and  also  toward  the  displaced  transmit  antenna.  As  a  result  the  Doppler 
spectrum  is  strongly  skewed  in  the  opposite  direction  to  that  of  Run  1.  Note 
SK  has  correspondingly  changed  sign  and  increased  in  magnitude  over  Run  1. 

The  standard  deviation  of  the  skew  has  also  increased  so  that  the  ratio 
of  average  skew  to  its  standard  deviation  remains  close  to  4.  More  data 
would  be  required  to  determine  whether  or  not  this  ratio  tends  to  remain 
constant  with  displacement  of  the  composite  antenna  pattern.  Nevertheless, 
these  data  suggest  that  the  Skew  estimator  is  as  good  or  better  than  the  Slope 
estimator  as  an  alignment  indicator,  but  with  the  advantage  that  absolute 
clock  setting  accuracy  is  not  a  requirement  with  Skew.  The  Skew  approach  does 
depend  on  the  extra  step  of  setting  the  clock  to  make  oi  equal  to  zero  before 
each  Skew  measure  is  attempted.  This  step  is  not  time  consuming  and  we 
estimate  from  this  experience  that  complete  Step  1  and  Step  2  alignment  pro¬ 
cedures  could  be  completed  in  about  ten  minutes  with  the  TRC-97  as  configured 
in  these  experiments.  Because  of  time  limitations  it  was  not  possible  to 
attempt  Step  2  of  the  alignment  procedure.  There  is  no  reason  to  believe  that 
the  results  would  differ  greatly  from  those  experienced  with  the  slope  algorithm 
since  ratios  of  mean  to  standard  deviation  of  about  4  were  observed  with  both 
Slope  and  Skew  measurements. 

Overview  of  Test  Results 

Both  the  "slope”  and  "skew"  estimation  techniques  were  shown  to  provide  a 
sensitive  measure  of  antenna  misalignment  and  a  fast  and  accurate  method  of  alignment 
based  on  these  estimates  was  demonstrated.  The  "correlation"  technique  could  not  be 
tested  because  the  program  was  devised  to  eliminate  fast  fading  data  as  a  safeguard 
against  aircraft  signals.  The  fadinq  rate  on  the  RADC  link  was  so  fast  that  the 
scheme  was  nearly  always  defeated.  However,  both  "slope"  and  "skew"  were  sufficiently 
successful  to  perhaps  negate  the  need  to  test  the  correlation  approach. 
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Accuracy  and  Alignment  Time 

Both  the  successful  modes  of  operation  provided  azimuth  alignment  to  within 
+_  1/3  degree  of  the  great  circle,  which  is  about  the  same  as  the  backlash  on  the 
TRC-97  antenna  gear  drive.  This  degree  of  accuracy  is  implied  by  the  fact  that 
at  1/3  degree  error,  the  average  value  of  the  displayed  quantity  equaled  the  rms  value 
of  its  1-sec  fluctuation.  The  displayed  quantity  is  derived  from  16  seconds  of  data, 
and  this  is  the  amount  of  time  required  for  the  rms  fluctuation  to  reach  their 
minimum  value  after  data  taking  is  started.  It  is  important  to  point  out  once  again, 
that  an  operator  of  the  aligner  would  tend  to  mentally  average  the  1-second 
displayed  numbers  before  making  an  estimate  of  the  misalignment.  If  15  seconds  of 
observations  are  permitted,  a  total  of  30  seconds  would  elapse  before  each 
adjustment  of  the  beam  could  be  made.  Assuming  about  ten  such  adjustments  are 
necessary  to  reach  final  alignment,  300  seconds  of  data  processing  time  would  be 
required  to  align  one  antenna  using  the  slope  technique.  Data  averaging  by  the 
operator  will  also  tend  to  improve  the  alignment  accuracy,  so  it  would  be  reasonable 
to  predict  that  accuracy  at  least  twice  as  good  as  +_  1/3°  could  be  achieved  in 
practice. 

Comparison  of  Slope  and  Skew  Methods 

The  tests  indicated  that  the  slope  and  the  skew  methods  achieved  essentially 
the  same  degree  of  accuracy.  However,  the  slope  method  is  not  usable  unless  some 
method  is  available  to  synchronize  the  transmitter  and  receiver  clocks.  If  the 
propagation  medium  provides  a  strong  specularly  reflected  signal  with  zero  Doppler 
shift  the  aligner's  clocks  can  be  synchronized  by  zeroing  the  average  frequency 
read  on  the  display  using  the  aligner's  "clock"  mode.  In  the  absense  of  a  specular 
signal  it  is  necessary  to  use  the  skew  mode  and  to  correct  the  receiver  clock  at 
each  new  antenna  position.  Correcting  the  clock  takes  an  additional  10  or  15 
seconds  at  each  antenna  position  which  about  doubles  the  time  required  to  align 
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the  antennas.  An  estimate  of  the  total  time  to  align  using  the  s.kew  mode  is  thus 
about  600  seconds. 

Our  experience  at  troposcatter  frequencies  near  1  GHz  indicates  specular  layers 
are  usually  present.  At  the  TRC-97  frequency  near  5  GHz  specular  reflections  are 
much  less  likely  to  occur.  Based  on  this  the  skew  algorithm  is  the  only  choice  for 
the  TRC-97  radios  providing  alignment  is  about  10  minutes  for  each  antenna.  Newer 
coherent  radios  such  as  the  TRC-170  with  atomic  frequency  control  would  be  able  to 
use  the  slope  algorithm  and  achieve  alignment  in  about  5  minutes  for  each  antenna. 

Connecting  the  Aligner  to  the  TRC-97 

The  overall  block  diagrams  in  Fig.  C.l  shows  the  aligner  attached  to  the  TRC-97 
transmitter  and  receiver  when  the  system  is  in  the  alignment  mode.  This  mode  is 
accomplished  by  placing  the  align-nonal ign  front  panel  switch  in  the  align 
position  which  applies  the  IF  output  to  the  aligner's  product  detectors  and 
substitutes  the  aligner's  coherent  frequency  sources  for  the  TRC-97 's  incoherent 
sources.  The  operating  frequency  in  the  alignment  mode  is  4880  MHz.  Figure  C.2a 
and  c.2b  show  details  of  the  relay  switching  of  the  signal  lines  from  the  TRC-97 
to  the  aligner  and  back.  In  Fig.  C.2a  the  port  labeled  "INCOHERENT  50  MHz  IN" 
receives  the  TRC-97's  50  MHz  synthesizer  output  that  normally  connects  to  J3A17. 

The  port  labeled  "INCOHERENT  70  MHz  IN"  receives  the  TRC-97 's  50  MHz  synthesizer 
output  that  normally  connects  to  Jl  of  the  upconverter  A17A4.  The  remaining  signal 
lines  and  ports  are  self-explanatory  and  should  be  followed  in  installing  the 
aligner  with  the  TRC-97. 

When  the  align-nonalign  switch  is  in  the  nonalign  position  or  when  the  aligner 
is  powered-off  the  TRC-97  is  returned  to  its  normal  configuration. 

Step-by-Step  Operating  Procedure  for  Skew  Method 

1.  Test  the  aligner's  microprocessor  by  depressing  the  "test"  pushbutton.  The 

-2 

display  should  read  019(-2)  which  represents  the  decimal  number  19  x  10  .  If 

it  reads  correctly  proceed  to  step  2. 
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2.  With  the  aligner  installed  according  to  the  preceding  description  and  the 

al ign-nonalign  switch  in  the  align  position  the  tropo  signal  should  be  acquired 
following  the  standard  procedures  outlined  in  the  TRC-97  instruction  manual.  If 
desired^  the  inphase  and  quadrature  outputs  from  the  back  panel  of  the  aligner 
can  be  displayed  on  an  oscilloscope  to  aid  in  visualizing  received  sional 
strength. 

1.  With  both  antennas  focused  the  aliqner  should  be  placed  in  the  clock  mode  by 
depressing  the  front  panel  pushbutton  labeled  "clock.”  The  displayed  average 
frequency,  to,  should  then  be  set  to  zero  with  the  front  panel  clock  frequency 
adjustment  control. 

4.  With  to  set  at  zero  the  skew-estimate  is  made  by  depressing  the  "skew"  pushbutton. 
The  numbers  displayed  will  vary  randomly  each  second  around  an  average  value. 

(If  an  aircraft  intrudes  the  common  volume  the  skew  measure  will  vary  widely 
with  time  as  the  aircraft  Doppler  varies.  The  "reset"  pushbutton  will  clear 
the  memory  allowing  a  new  skew  estimate  to  be  obtained  when  the  aircraft  clears 
the  antenna  beams.) 

5.  When  a  stable  skew  estimate  is  obtained  the  receive  antenna  should  be  rotated 
about  1/2°  in  one  direction  and  fixed. 

6.  The  clock  mode  is  again  used  to  set  (0  equal  to  zero. 

7.  The  skew  mode  is  next  used  to  reestimate  the  skew.  If  the  average  value  has 
been  reduced  the  correct  direction  of  antenna  rotation  was  chosen.  Continue  to 
reset  the  antenna  in  small  increments  until  the  skew  display 

changes  sign.  This  is  evidence  that  the  receiving  beam  has  gone  through  the 
symmetrical  stradding  beam  position. 

8.  With  the  receiving  beam  positioned  to  give  zero  average  skew,  the  defocuser  is 
attached.  If  the  skew  remains  zero  the  beams  are  correctly  aligned.  If  the 
skew  becomes  non-zero  proceed  to  step  9. 
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9.  The  transmitter  and  receiver  beams  are  now  swung  synchronously  in,  say,  1/4° 
increments,  both  clockwise  or  both  counter  clockwise.  An  order-wire  voice 
channel  back  to  the  transmitter  is  necessary  to  coordinate  this  synchronous 
beam  swinging  procedure.  At  each  position  U)  is  set  to  zero  and  the  skew 
measured.  If  the  correct  direction  of  rotation  (CW  or  CCW)  has  been  chosen, 
the  skew  measure  will  fall  with  each  new  increment  until  it  reaches  zero.  At 
this  point  the  two  beams  are  producing  a  non-skewed  Doppler  spectrum  indicating 
that  they  have  met  at  the  great-circle. 

10.  The  receiving  beam  is  then  refocused  to  complete  the  alignment  procedure. 

11.  It  should  be  noted  that  if  the  incorrect  direction  of  rotation  had  been  chosen 
the  skew  would  have  increased  as  the  beams  diverged  from  the  great-circle. 

Step-by-Step  Operating  Procedure  for  the  Slope  Method 

1.  The  signal  should  be  acquired  by  following  the  standard  methods  outlined  in  the 
TRC  instruction  manual. 

2.  The  atomically  controlled  5  MHz  receiver  clock  should  be  substituted  for  the 

5  MHz  quartz  standard  contained  in  the  aligner.  It  is  also  assumed  that  the 

If 

transmitter  and  receiver  atomic  standards  are  in  zero-beat. 

3.  The  aligner's  microprocessor  is  tested  by  depressing  the  "test"  pushbutton. 

-2 

The  correct  display  is  019-2  or  19  x  10 

4.  Depressing  the  "slope”  pushbutton  will  cause  the  aligner  to  read  the  slope 
of  the  crosscorrelation  function  of  the  quadrature  components  of  the  received 
siqnal.  The  slope  will  be  non-zero  when  the  antennas  are  misaligned. 

5.  The  receiving  beam  should  now  be  moved  in  azimuth  about  1/2  degree  and  fixed. 

6.  The  slope  is  remeasured.  If  it  has  decreased  the  antenna  movement  was  in 
the  correct  direction. 

7.  Continue  to  increment  the  antenna  in  the  direction  that  decreased  slope  until 
the  slope  equals  zero. 
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H.  The  receiving  antenna  is  next  dofocused.  If  the  slope  remains  zero  both  beams 
are  correctly  aligned.  If  the  slope  becomes  non-zero  proceed  to  step  9. 

9.  The  transmitting  and  receiving  beams  are  now  swung  synchronously,  both  clockwise 
or  both  counter-clockwise  in  say,  1/4  degree  steps.  If  the  correct  direction 
is  chosen  the  slope  will  decrease  with  each  increment  until  it  reaches  zero. 

At  this  point  the  beams  are  correctly  aligned. 

10.  The  defocuser  is  removed  to  complete  the  alignment  procedure. 

Suggestions  For  Further  Study 
1 .  Deficiency  of  Skew  Estimation  Program 

For  reasons  that  are  not  entirely  clear,  the  skew  algorithm  displays 
sensitivity  to  clock  error.  As  noted  in  the  introduction  this  is  contrary 
to  theoretical  expectations,  but  may  be  due  to  the  fact  that  the  Doppler 
spectrum  is  bounded  by  an  8  Hz  low-pass  digital  filter.  Therefore  when 
the  X  (t)  or  Y  ft)  spectrum  is  spread  by  offsetting  the  clock,  either  the 
upper  or  lower  edges  of  the  Doppler  spectrum  are  attenuated  depending  on 
whether  the  clock  is  raised  or  lowered  in  frequency  from  zero  beat.  This 
could  induce  positive  or  negative  skew  respectively,  widening  the  filter 
bandwidth  should  be  tried  as  a  solution  to  this  problem. 

Another  factor  that  should  be  considered  is  the  possibility  that 
the  skew  programming  of  the  microprocessor  is  deficient  in  computing  the 
first  and  second  derivatives  of  the  quadrature  components  of  the  received 
signal.  Since  these  signals  are  first  quantized  by  an  A/D  converter, 
successive  samples  falling  in  the  same  level  bin  will  yield  a  first 
derivative  of  zero  and  three  successive  samples  in  the  same  bin  will  yield 
a  second  derivative  of  zero.  Thus,  when  the  signal  derivatives  tend  to  be 
small,  they  may  often  be  under  estimated.  Slowly  fluctuating  signals  may 
therefore  lead  to  errors  in  the  various  terms  that  make  up  the  Skew 
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equation.  This  equation  yields  a  small  difference  from  large  numbers  and 
is  therefore  quite  sensitive  to  errors  associated  with  the  various  terms. 

See  equation  (19). 

During  UW  tests  Doppler  spreads  were  much  narrower  than  on  the  RADC 
link.  We  experienced  skew  estimates  Sensitive  to  clock-offset  particularly 
when  the  Doppler  spectrum  was  extremely  narrow  as  evidenced  by  a  zero-Doppler 
spike  in  the  spectrum,  suggesting  the  derivatives  as  the  primary  cause.  The 
8  Hz  filter  cutoff  apparently  was  no  problem  on  the  UW  link,  since  we  were 
able  to  obtain  skew  estimates  insensitive  to  clock-offset  on  two  occasions 
when  the  Doppler  spectrum  was  narrower  than  8  Hz,  but  free  from  a  specular 
component . 

If  the  above  two  problems  prove  incorrect  as  the  cause  of  skew  sensitivity 
to  clock  offset,  it  might  be  expedient  to  attempt  to  make  clock  correction  (to 
force  0)  to  equal  zero)  an  automatic  process.  Since  to  is  outputted  by  the 
slope  algorithm  an  estimate  of  it  could  be  used  to  correct  the  receiver  clock 
by  feedback  to  the  clock  control  voltage  input.  This  would  remove  the  extra 
step  that  otherwise  must  be  provided  by  the  operator. 

2 .  Additional  Testing  Requirements 

Time  permitted  a  fairly  limited  amount  of  testing  on  the  RADC  link. 

During  these  tests  the  fading  rate  was  high  and  signal  characteristics  were 
consistent  for  the  several  days  of  the  tests.  Testing  under  a  variety  of 
signal  conditions  should  be  attempted  before  decisions  concerning  filter 
bandwidths  are  finalized. 

The  skew  algorithm  was  tested  only  with  focused  beams.  It  is  reasonable 
to  expect  the  final  accuracy  with  a  defocused  beam  to  parallel  the  results 
of  the  slope  algorithm  test,  however.  Nevertheless,  it  would  be  appropriate 
to  check  this  when  additional  tests  are  run  under  different  signal  conditions. 
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Fig.Al  BLOCK  DIAGRAM  OF  MICROPROCESSOR  SYSTEM 
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A.  14 


Parts  List  Microprocessor  Board  #1 


Location 


1  tem 


Mil 
KH 
14 
[•'2 
C'2 
K16 
1117 
K17 
C26 
K  37 
H  39 
C3R 
A3R 
M4R 
K4R 
148 

f  JC,  ■> 

C49 
A49 
E  39 
K65 
M64 
AS  9 
1161 
IP.  1 
CIS 
C15 
CIS 
CIS 
CIS 
1161 
H61 
F.61 
E61 
K16 
M27 


Intel  2114,  1C 
Intel  2114,  IC 
Intel  2114,  1C 
Intel  2716,  Ic 
Intel  8212,  Ic 

CTS  760-1-R1K,  Resistor  Network 
Intel  81SS,  Ic 
Intel  808S 

National  DM  74I.il 32,  IC 

Nat  ional  I'M  74l.it  1  3R,  IC 

National  DM  741.POR,  IC 

National  DM  74I.R13R,  IC 

National  DM  74I.S155,  IC 

National  DM  74I.S244,  1C 

National  DM  74LS244,  IC 

National  DM  74LS273,  IC 

National  DM  741.S74,  IC 

National  DM  74LS04,  1C 

National  DM  7406,  IC 

American  Micro  AM9S11,  IC 

National  AM  3705CN,  IC 

Analog  Devices  ADS71JD,  \/D  converter 

Intel  2114,  IC 

National  LF39R,  S/ll  IC 

National  LF398,  S/H  IC 

56K  Resistor  1/4  W,  S% 

4  MHz  Crystal 

180  pF  Capacitor-Ceramic 

1  li F  Capacitor 

1N44  Diode 

.01  ;i F  Capacitor-Ceramic 
0.1  liF  Capacitor-Ceramic 
.01  li F  Capacitor-Ceramic 
0.1  li F  Capacitor-Ceramic 
IK  Resistor  1/4  W,  5% 

Dip  Switch  8  section,  SPST 


Tarts  List  Microprocessor  Board  #2 


KS 

Intel 

2114, 

IC 

15 

Intel 

2114, 

IC 

05 

Intel 

2114, 

IC' 

E5 

‘Intel 

2114, 

IC 

717 

Intol 

2716, 

IC 

Cl  7 

Intel 

8253, 

IC 

El  7 

Motorola  MC14R9,  IC 

E27 

Motorola  MC 

1488, 

IC 

7  32 

Into  1 

8251 , 

IC 

0  32 

Intel 

8251 , 

IC 

E  37 

National  DM74LS14, 

IC 

('37 

National  DM74I.S01, 

IC 

047 

National  DM74I.S32, 

IC 

749 

Intel 

2716, 

IC 

N42 

IK  Resistor 

1/4  W, 

5% 

M49 

4.7  K 

Resistor  1/4 

W, 

N49 

4.7K 

Resistor  1/4 

W, 

Tarts  List  Microprocessor  Display 

1  -  National  DM74LS156,  IC 

1  -  National  DM74LS175,  IC 

7  -  Monsanto  MAN72,  Display 

7  -  Fairchild  9374,  IC 


Power  Supplies 


+15  Volt,  100  mA  Tower  Mato  MD15A 
+  5  Volt,  6  A  Power  Mate  ETR-142E 


APPENDIX  B  -  MICROPROCESSOR  PROGRAMMING 
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APPENDIX  C 


SYNTHESIZER  DRAWINGS  AND  PARTS  LISTS 


SYNTHESIZER 


70 

KHi 


DETECTORS 


ALIGNER  CONNECTIONS  TO  RECEIVER 


J 


ALIGNER  CONNECTIONS  TO  TRANSMITTER. 


J 


Figure  C.l  INTERCONNECTIONS  OF  ALIGNER  AND  TRC-97 


HoiOBina  Lb 


SIGNAL  LEAD  RELAY  SWITCHING  FOR  TRC-97  RECEIVER  APPLICATION 


5  MHz 

FREQUENCY 

STANDARD 


Figure  C.4. 


TRANSMITTER  FREQUENCY  SYNTHESIZER  BASIC  DESIGN 


DIAGRAM 


gure  C.10.  50  MHz  AMPLIFIER 


o 


Figure  C.llb.  LOW  PASS  AND  BANDPASS  FILTER 


VOLTS  REGULATED  TO 


NEGATIVE  VOLTAGE 


Figure  C.13 


50  MHz  Amplifier  Output  Spectrum.  Full  Scale  =  20  dbm; 
Vertical  Scale:  10  db/Div;  Horizontal  Scale:  10  MHz/Di v 
Center  Freq.  =  0  MHz. 


t 


Figure  C.14.  70  rillz  Amplifier  Output  Spectrum.  Full  Scale  =  20  dbm; 

Vertical  Scale:  10  db/Div;  Horizontal  Scale:  10  MHz/Di 
Center  Freq.  =  0  MHz. 
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Figure  C.I5b.  PHASE  LOCK  LOOP  GROUND  PLANE 
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70  MHz  PLL 


CI00.  101.  108.  116 
L 100.  L 102.  R 100 


\ 


ON  WALL  OF  BOX 


Figure  C-19.  70  MH-s  Phase  Lock  Loop.  The  50  MHz  PLL  is  identical 

except  for  part  numbers  in  the  two-hundreds. 


70  MHz  AMP 


C3I0.  309  1 

R  307  I  MOUNTED  ON  WALL 


Figure  C-20.  70  MHz  amplifier. 


50  MHz  AMP 

C4 1 1.  R407,  C4  ON  BOX  WALL 


C5 13.  512,  SOS.  RSIO.  RSI  I 
MOUNTED  ON  WALL  OF  BOX 


Figure 


C-22.  Detector  and  Low-Pass  Filter 
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Parts  List  70  MHz  PLL 


Capacitors : 

C100  1200  pF 

C101  1200  pF  Feedthrough 

C102  10  microFarad,  25  Volt, Electrolytic 

C103  0.01  microFarad 

C104  330  pF 

C105  120  pF 

C106  120  pF 

C107  selected 

C108  1200  Feedthrough 

C109  0.01  microFarad 

C110  330  pF 

Clll  330  pF 

C112  330  pF 

C113  330  pF 

C114  20  pF 

C115  330  pF 

C116  1200  pF 

C117  47  microFarad,  25  Volt,  Electrolytic 

C118  330  pF 

C119  330  pF 

C120  330  pF 

C121  330  pF 

C122  330  pF 

C123  4.7  pF 

Transistors : 

Q100  2N918 

Q101  2N3866 

Resistors:  All  \  Watt,  5% 

R100  500  Ohm,  20  turn  potentiometer 

R101  470  Ohm 

R102  2k  Ohm 

R103  2k  Ohm 

R104  2k  Ohm 

R10 5  5.1k  Ohm 

RIOS  5.1k  Ohm 

R107  510  Ohm 

R108  20k  Ohm 

R109  33k  Ohm 

R110  18k  Ohm 

Rill  2.2k  Ohm 

R112  2.7k  Ohm 

R113  lk  Ohm 

R114  lk  Ohm 

R115  12k  Ohm 

R116  47k  Ohm 

R117  330  Ohm 


Diodes : 

D100  1N4148 

D101  1N4148 
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Inductors : 
L100  4.7 
L101  4.7 
L102  4.7 


microHenry 

microHenry 

microHenry 


Integrated  Circuits: 

U100  XR215 

U101  74S132 

U102  93S05 

U103  Mini  Circuits  Laboratory  SRA-1  balanced  mixer 


B 


Capacitors : 

C200  1200  pF 

C201  1200  pF  Feedthrough 

C202  10  microFarad,  25  Volt,  Electrolytic 

C203  0.01  microFarad 

C204  330  pF 

C205  120  pF 

C206  120  pF 

C207  selected 

C208  1200  pF  Feedthrough 

C209  0.01  microFarad 

C210  330  pF 

C211  330  pF 

C212  330  pF 

C213  330  pF 

C214  43  pF 

C215  330pF 

C216  1200  pF 

C217  47  microfarad,  25  Volt,  Electrolytic 

C218  330  pF 

C219  330  pF 

C220  330  pF 

C221  330  pF 

C222  330  pF 

Transistors : 

Q200  2N918 

Q201  2N918 

Resistors:  All  %  Watt,  5^ 

R200  500  ohm,  20  turn  potentiometer 

R201  470  Ohm 

R202  2k  ohm 

R203  2k  ohm 

R204  2k  Ohm 

R205  5.1k  Ohm 

R206  5.1k  Ohm 

R207  510  Ohm 

R208  20k  Ohm 

R209  33k  Ohm 

R210  18k  Ohm 

R211  2.2k  Ohm 

R212  2.7k  Ohm 

R213  lk  Ohm 

R214  lk  Ohm 

R215  12k  Ohm 

R216  47k  Ohm 

R217  330k  Ohm 

Diodes : 

D200  1N4148 

D201  1N4148 
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Inductors : 

L200  4.7  microHenry 

L201  4.7  microHenry 

L202  4.7  microHenry 

Integrated  Circuits: 

U200  XR215 

U201  74S132 

U202  93S05 

U203  Mini  Circuits  Laboratory  SRA-1  balanced  mixer 

U204  Mini  Circuits  Laboratory  PSC-2  power  splitter 


Parts  List  70  MHz  Amplifier 


Capacitors : 
C300  1000  pF 

C301  62  pF 

C302  1000  pF 

C303  16  pF 


C30  5  1000  pF 

C307  18  pF 

C308  1000  pF 

C309  1200  pF  Feedthrough 

C310  1200  pF 

C311  47  microFarad,  25  Volt,  Electrolytic 

C312  0.01  microFarad 

C314  12  pF 

C315  1000  pF 

C316  1000  pF 

Resistors:  All  \  Watt,  5% 

R300  100  Ohm 

R301  2.2k  Ohm 

R302  3.9k  Ohm 

R303  330  Ohm 

R304  2.2k  Ohm 

R305  2.7k  Ohm 

R306  330  Ohm 

R307  10  Ohm 

Inductors : 

L300  0.28  microhenry 

L301  2.2  microHenry 

L302  2.2  microHenry 

L303  0.28  microHenry 

Transistors : 

Q300  2N918 

Q301  2N918 


Parts  List  50  MHz  Amplifier 


C. 


Capacitors : 

C400  27  pF 

C401  62  pF 

C402  1000  pF 

C403  1000  pF 

C404  1000  pF 

C405  27  pF 

C406  47  microFarad,  25  Volt,  Electrolytic 

C407  0.01  microFarad 

C408  1000  pF 

C409  1000  pF 

C410  1000  pF 

C41 1  1200  pF  Feedthrough 

C412  1200  pF  Feedthrough 

Resistors:  All  h  Watt,  5^ 

R400  100  Ohm 

R401  2.2k  Ohm 

R402  2.7k  Ohm 

R403  330  Ohm 

R404  3.9k  Ohm 

R405  2.2k  Ohm 

R406  120  Ohm 

R407  10  Ohm 

Inductors: 

L400  0.39  microHenry 

L401  0.39  microHenry 

L402  2.2  microHenry 

L403  0.39  microHenry 

L404  2.2  microHenry 

Transistors : 

Q400  2N918 

Q401  2N3866 
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Parts  List  Detector 


Capacitors : 

C500 

0.001  microFarad 

C501 

0.1  microFarad 

C502 

0.1  microFarad 

C503 

0.001  microFarad 

C504 

0.1  microFarad 

C505 

0.1  microFarad 

C506 

0.01  microFarad 

C507 

100  microFarad,  15  ' 

C508 

1200  pF  Feedthrough 

C509 

1200  pF 

C510 

0.01  microFarad 

C511 

100  microFarad 

C512 

1200  pF  Feedthrough 

C513 

1200  pF 

C514 

0.01  microFarad 

C515 

0.01  microFarad 

C.'!>16  0.01  microFarad 

Cr>17  0.01  microFarad 

Resistors:  All  Watt,  5 % 

R500 

51  Ohm 

R501 

820  W 

R502 

i.’K 

R503 

10k  Ohm 

R504 

100k  Ohm 

R505 

51  Ohm 

R506 

ik  ■; 

R507 

i.'k  s: 

R508 

10k  Ohm 

R509 

100k  Ohm 

R510 

R511 

R512 


10  Ohm 
10  Ohm 
51  Ohm 


Integrated  Circuits: 

U500  Mini  Circuits  Laboratory  PSCQ2quadrature  power  splitter 

U501  "  "  •»  SRA-1  balanced  mixer 

U502  "  "  "  SRA-1  balanced  mixer 

U503  "  "  ••  PSC2-1  power  splitter 

U504  LM741 

U505  LM741 


Operational  Amplifier;;  U601-U610  13741  (LM  741) 

U611  LM  7815  +15  volt  regulator 
U612  LM  7915  -15  volt  regulator 


Capacitors  C619,  C620  100  microfarad,  56  volt  DC 

Capacitors  C601  -  C618,  1  microfarad,  50  volt 


Resistors  $641,  642,  150  12,  10  watt 


All  resistors  arc 

1/4  watts 

R601,  R614, 

R627 , 

R634 

— 

313 

kfi 

R602 ,  R605 , 
R615,  R616 , 

R606 , 
R619, 

R707  'l 
R620  j 

— 

510 

k!2 

R60  3 ,  R617 

— 

12 

k!2 

R628 ,  R6  35 

— 

15 

k£2 

R604 ,  R609, 

R613 , 

R618  > 

R622 ,  R626 , 
R636 ,  R639 

R629, 

R633 

10 

kfl 

R608 ,  R621 

— 

1.5 

kQ 

R610 ,  R611, 
R630,  R6  31 , 

R62  3 , 
R6  37 , 

R624  \ 
R638  J 

— 

6.8 

k!2 

R612 ,  R625 , 

R632  , 

R640 

— 

5.6 

k!2 

t 
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Capacitors : 

C vOO  1000  microFarad,  56  Volt,  Electrolytic 

•C701  1000  microFarad,  56  Volt,  Electrolytic 

C  702  1200  pF  Feedthrough 

Resistors : 

R  700  100  Ohm,  10  Watt 

R 7oi  150  Ohm,  10  Watt 

•R  702  2  00  Ohm,  5  Watt 

Rectifiers : 

CR 700  2  Ampere,  50  Volt  PIV,  Diode  Bridge 
•CR701  2  Ampere,  50  Volt  PIV,  Diode  Bridge 

Integrated  Circuits: 

U7D0  LM7824,  24  Volt  Regulator 

U701  LM7812,  12  Volt  Regulator 

U  702  LM7805,  5  Volt  Regulator 

•U703  LM7912,  -12  Volt  Regulator 

T700  Power  Transformer,  Primary  Voltage  =  117  VRMS,  60  Hz 

Two  secondaries,  28  VRMS,  1  Ampere  each 
eg.  Signal  Transformer  Co.  Inc.  Model  56-1 

•Not  required  for  the  transmitter  power  supply. 
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Interrelations  of  Synthesizer  Drawings  and  Parts  Lists 


Module 

Schematic 

Printed  Board 

Parts  Layout 

Parts  List 

70  MHz  PPL 

Fig . 

C.  7 

Fig.  C.15a,b 

Fig.  C . 19 

Page  C29 

50  MHz  PPL 

Fig. 

C .  8 

Fig.  C.15a,b 

Fig.  C. 19 

Page  C31 

70  MHz  Amp 

Fig. 

C.  9 

Fig.  C.17a,b 

Fig.  C. 20 

Page  C33 

50  MHz  Amp 

Fig. 

C.10 

Fig.  C.10a,b 

Fig.  C.21 

Page  C34 

Dotector-Fi Iter 

Fig. 

C.  11a 

Fig.  C.18a,b 

Fig.  C . 22 

Page  C35 

LP-BP  Filter 

Fig. 

C.llb 

Fig.  not  used 

Fig.  C. 23 

Page  C36 

Power  Supply 

Fig. 

C .  1 2 

Fig.  not  used 

Page  C37 

50  MHz  Spectrum 

Fig. 

C.  13 

70  MHz  Spectrum 

Fig. 
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Appendix  D 


Theory  of  Doppler  Alignment  Estimators 
and 

Tropo-Channel  Simulation 
Contents 


Tart  I 

"Doppler  Spectra  as  a  Function  of  Tro[>o-Link  Parameters  and 
Accuracy  Analysis  of  'Slope'  Estimators" 

Interim  Report  for  12/11/77  to  1/11/78 . .....D,2 

Part  II 

"Estimation  of  Quadrature  Correlation  Function  for  all  Delays  by 
Analog  and  Polarity-Coincidence  Methods” 

Interim  Report  for  1/11/78  to  2/11/78 . . . D.26 

Part  III 

"A  New  Complex  Gaussian  Troposcatter  Channel  Simulator" 

Interim  Report  for  6/11/78  to  7/11/78 . D.49 

Summary  of  Appendix  D 

Part  I  reviews  the  relationship  between  troposcatter  link  geometry,  the 
meteorology  of  the  scattering  structure,  and  the  Doppler  spectrum  imparted  by 
the  mean  cross-path  wind.  It  is  shown  that  Doppler  spectrum  symmetry  or  its 
lack  can  be  used  to  test  for  correct  antenna  alignment.  Various  estimators 
for  measures  of  symmetry  are  then  discussed  including  average  Doppler  shift 
and  skew.  The  variance  of  these  estimators  and  their  sensitivity  to  link 
parameters  are  calculated  and  tabulated  for  several  types  of  signal  processing. 

Part  II  is  essentially  devoted  to  the  analysis  of  the  variance  of  the 
cross-correlation  or  slope  estimator  when  one  or  both  quadrature  component 
signals  are  hard  clipped.  The  increase  in  variance  of  6  db  due  to  clipping 
was  the  main  reason  for  using  non-distorted  signals  in  the  final  version 
of  the  aligner. 

Part  III  is  a  description  of  a  tropochannel  simulator  in  which  delay-line 
tap  gains  are  Gaussian  random  variables  with  quadrature  statistics 
that  produce  Doppler  effects  analogous  to  those  found  in  true  channels. 

Skewed  spectra  can  be  synthesized  and  thus  the  effects  of  aircraft  intrusion 
can  be  modeled. 


]). 
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Doppler  Spectra  as  a  function  of  Tropn-link  Parameters 
and  A<  apacy  Analysis  of  ‘Slope1  Estimators 


Introduction 

When  a  monochromatic  carrier  is  transmitted,  the  received  signal  is 
normally  spread  about  the  carrier  frequency  because1  of  atmospheri cal ly- 
induced  Doppler  shifts.  In  this  chapter,  the  relative  accuracy  of  the 
alignment  scheme  by  'slope  method'  is  related  to  the  Doppler  spread  function 
which  is  shown  to  depend  on  the  antenna  patterns  and  their  position  relative 
to  the  great-ci rcle ,  the  angular  dependence  of  scattering,  the  cross  path 
wind  speed,  and  the  link  geometry.  Next  we  show  that  the  accuracy  of  the 
slope  estimators  is  degraded  by  factors  which  increase  the  fade-rate  such 
as  beam-width  or  cross-path  wind  speed.  Finally,  we  show  that  the 
accuracy  is  further  degraded  if  'hard-clipped'  estimator  is  used  instead 
of  'true'  estimator  of  slope  of  quadrature  correlation  function  at  zero  lag. 

Review  of  Link  Geometry:  The  paths  taken  by  the  transmitted  and  scattered 
signals  are  slightly  bent  toward  the  earth  due  to  the  average  decrease  in 
refractive  index  with  height.  Bean  and  Dutton  (1966)  have  shown  that  under 
normal  conditions  the  average  refractivity  variation  near  the  earth's 
surface  is  linear,  thereby  permitting  the  use  of  an  earth  model  whose 
radius  is  approximately  4/3  that  of  the  true  earth  radius.  The  equivalent 
earth  model  has  several  distinct  practical  advantages:  1)  the  refractive 
index  surrounding  the  earth  is  uniform  and  constant,  2)  the  curved  radio 
paths  become  straight  lines  (or  ray-paths),  and  3)  the  ray  paths  have  the 
same  length  as  their  curved  counterparts.  Assuming  this  equivalent  earth 
model,  the  Doppler  spectrum  is  written  in  terms  of  angular  dependence  and 
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antenna  function  using  the  Doppler-shi ft  model  of  Birkemeier  (1968). 

To  each  point  in  space  there  corresponds  a  ray  path  length  L(x,y,z) 
which  is  defined  by  a  ray  from  the  transmitter  to  the  receiver  passing 
through  the  given  point.  Surfaces  of  constant  path  or  phase  are  approximately 
ellipsoids  of  revolution  having  the  radio  terminals  as  foci  (Fig.  2).  A 
Signal  obtained  via  any  single  moving  point  is  shifted  in  frequency  by  an 
amount  proportional  to  the  rate  at  which  the  point  cuts  equiphase  surfaces. 
Thus  the  'single  scatterer'  Doppler-shi ft  f  is  given  by 


f 


1  dL 
x  dt 


(1) 


where  x  is  the  wavelength  and  v  is  the  velocity  of  the  moving  scatterer. 
Using  the  coordinate  system  of  Fig.  1  with  origin  at  the  midpoint  of  the 
straight  (chordal)  line  between  the  terminals: 

L  =  [(d+x)2+y2+z2]%  +  [(d-x)2+y2+z2]%  (2) 

Applying  (1)  and  (2)  yields,  Birkemeier  (1968), 


f  = 


4  ux(l-4d2/L2)+vy+wz 
^  l-16x2d2/L4 


(3) 


where  (u,v,w)  are  the  (x,y,z)  components  of  v .  The  relative  importance 
Of  the  various  terms  in  (3)  depends  upon  the  position  of  the  scatterer 
and  on  the  relative  magnitude  of  the  velocity  components.  The  along  path 
(u)  contribution  is  negligible  because  the  ellipsoids  are  longitudinally 
vary  'flat'.  The  factor  (1  -  ^-)  is  small,  and  equal  to  4  z^2/!.2  =  4  hm2/d2 
where  z[n  and  1^  are  the  midpath  heights  of  an  ellipsoid  above  the  chordal 
axis  and  spherical  earth,  respectively.  This  factor  ranges  fru m  10*4  to 
_2 

10  from  the  bottom  to  the  top  of  the  common  volume.  Since  at  a  given 
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point,  the  average  horizontal  velocity  normally  greatly  exceeds  the 
average  vertical  velocity,  Eq.  (3)  can  be  simplified  to 

f  =  -  I  v(^)  =  -  j  v  sina  =  -  y  va  (4) 

where  a  is  the  azimuthal  angle  of  the  single  scatterer  from  the  great 
circle.  For  a  given  crosspath  wind  speed,  pointing  the  antennas  to 
azimuthal  angles  away  from  the  great  circle  tends  to  favor  scatterers  at 
larger  angles  a,  which  consequently  increases  the  observed  Doppler  shift. 

There  are  many  scatterers  at  various  positions,  however,  that  create  a 
spectrum  of  Doppler  shifts  in  the  received  signal.  Birkemeier  (1968) 
specified  the  dependence  of  Doppler  spectrum  upon  azimuthal  pointing 
angle  a  in  terms  of  the  multiplied  azimuthal  pattern  function  G(a)  of 

a 

the  antennas  and  the  scattered  power  per  unit  azimuthal  angle  A(a).  For 
synchronously  pointed  antennas  with  boresights  at  a,, 

a 

S(f)  =  S  (a.aj  «  G(a-aJA(o)  (5) 

a  a  a 

where  S ( f )  is  the  Doppler  spectrum,  S  (a, a  )  is  Doppler  spectrum  written 

a  a 

in  terms  of  a.  And  a  and  f  are  related  by  Eq.  (4). 

Composite  Antenna  Function 

The  antenna  functions  used  in  the  simulation  are  contrived  to  be  a 
simple  but  accurate  representation  of  the  actual  antennas  used  in  the 
radio  experiments.  The  parameters  of  link  antennas  at  the  UW  link  between 
Arlington,  Wisconsin  and  Cedar  Rapids,  Iowa  and  RADC  link  between  Youngstown 
and  Verona,  New  York  are  given  in  Table  1.  Under  normal  operating  conditions, 
the  antennas  are  bore-sited  at  the  horizon.  Thus,  rays  from  the  lower  half 
of  the  beam  are  directed  at  the  foreground  and  are  both  absorbed  by  the 


terrain  and  diffusely  reflected.  Rays  from  the  upperhalf  of  each  beam 
are  directed  over  the  horizon  into  the  troposphere  and  their  intersections 
are  referred  to  as  the  common  volume. 

Representation  of  the  antenna  field  pattern  by  Gaussian  functions  is 
analytically  convenient  and  reasonably  accurate.  If  the  transmitting 
antenna  is  assumed  pointed  at  the  horizon  over  the  great  circle  but  the 
receiver  antenna  is  azimuthally  pointed  at  an  angle  a  (azimuthal  height  y  ) 
with  respect  to  great  circle  then  the  nominal  voltage  gain  at  the  receiver 
from  any  scattering  point  at  a  given  height  zQ  is  given  by 


2 

GR(y)GT(y)  «  ^  exp['  exp[-(z0-zgr)2/2a/] 


i  (y-yJ  o  o 

1  exp[-  —y-]  exp[-(z0-zgr)2/2aR2] 
R  ZoR 


(6) 


where  rR-rT  =  distance  to  scattering  point  from  receiving  and  transmitting 
antennas  respectively. 

2  2  . 
oR  .a-j.  =  variance  of  Gaussian  approximation  determined  by  antenna 

beam  width. 

Within  the  common  volume 

zgr  =  grazing  ray  height  above  link  axis,  and 


2  _  trRa3)2 
'  I41n2l 


(rTa3) 

(41n2) 


and  y  =  da 

a  a 


where  is  -3db  beam-width  of  a  single  antenna  pattern  and  2d  is  the  length 
of  the  link  axis  as  shown  in  Fig.  1. 


1).  (> 


Assuming  =  d,  at  midpath  and  grouping  the  constants  in 


Eg.  (6)  by  G0(z0>d),  antenna  function  G(a)  can  be  written  as 


G(a)  =  G(y)  =  [Gs(y)GT(y)]2=  Gn(zn.d)  expt-^]  exp[-  3 


O'  o ’ 


2 

2—3 


(7) 


where 


2  _ 


and  y  =  da  with  y  =  da 


0  =  T4TFT2T  --  -  -  Ja  -a 

Substituting  Eq.  (8)  into  Eq.  (7),  one  obtains: 


(8) 


2  q  )<■ 

G(a)  =  Gg(z0>d)  exp[-  exp[-  - ] 


19) 


2  2 

where  =  03/4102. 

In  the  above  discussion,  single  Gaussian  function  has  been  used  to 
represent  the  main  lobe  of  the  composite.  For  better  refinement  with 
side  lobe  structure,  bi-Gaussian  or  penta-Gaussian  function  (Merrill,  1968) 
can  be  used. 

Angular  Dependence 


A  major  ingredient  of  turbulence  scattering  theory  is  angular 
dependence.  The  signal  received  from  a  scatterer  is  a  function  of  the 
angle  between  transmitted  rays  incident  upon  the  scatterer  and  the  rays 
from  the  scatterer  to  the  receiver.  The  angle  involved  is  called  the 
"scattering  angle".  The  reflection  coefficient  of  the  layer  theories 
(Friis,  Crawford  and  Hogg,  1957;  Ducastel,  et  al.,  1960)  involves  the 
same  "reflection  angle".  Merrill  (1968)  called  it  "deflection  angle"  which 
does  not  explicitly  indicate  a  particular  mechanism. 


A  combination  of  turbulence  theory  and  wave  theory  results  in  the 
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scattering  cross  section 

o  =  -—  sin2  6  sin*p(e/2)  (10) 

x"p 

where  e  is  the  scattering  angle  and  8  is  the  angle  between  the  electric 
field  vector  and  the  ray  to  the  receiver,  e  is  also  a  function  of 
position  but  is  always  nearly  90°.  The  factor  C  includes  the  effect  of 
wavelength,  scale  of  turbulence  and  mean-square  refractivity  fluctuation 
intensity.  For  a  thin  layer  with  homogeneous  statistical  properties,  C  is 
assumed  to  be  independent  of  position.  The  power  law  p  depends  upon  the 
particular  version  of  turbulence  theory  advocated.  The  prominent  Kolmogorov- 
Tatarski  theory  (Tatarski ,  1961)  gives  p  =  11/3. 

Based  upon  the  scattering  cross  section  function  of  Eq.  (10),  one 
obtains  the  angular  dependence  function  A(a) 

At)  ■  A(y)  =  Ao(sin  e/2)~11/3  (11) 


where 


sin  e/2 


2  2 


at  midpath  (x=0)  and  Ao  is  a  constant 


depending  upon  the  wavelength  and  intrinsic  strenath  of  the  scatterers.  y  and 
a  are  related  as  in  Eq.  (8).  Equation  (11)  can  be  rewritten  as 


A(a)  =  Ao(a2+z2/d2) 


■11/6 


(12) 


Frequency  dependence  of  A(a)  can  be  obtained  from  the  relation  between 
a  and  f  as  in  Eq.  (4). 

1 


Approximation  of  Angular  Dependence 
for  Small  Antenna  Pointing  Angles 

For  small  antenna  pointing  angles,  the  angular  dependence  function 
can  be  approximated  by  the  Gaussian  function  of  the  form 

,  „  2  (y^+:0) 

A(a)  =  A(y)  =  CQ  exp[-  i  *  ^  ]  (1 

X  d 

where  CQ  and  tc  are  constants  depending  upon  the  3-dimensional  space 
correlation  length.  For  a  given  thin  layer  at  a  height  zQ,  one  can 
rewrite  Eq.  (13)  as 

2 

A(a)  =  A0(z0)exp[-  fj]  (1 

A 

2 

where  specifies  the  severity  of  rolloff  of  A(a)  and  is  dependent 
upon  wavelength  X. 


Functional  Representation  of  Doppler  Spectrum 
When  no  Clock-offset  is  Present 

From  Eqs.  (5),  (9)  and  (14)  one  can  write  the  spectral  density  of  the 
received  signal  as: 


2  (a_a  j  ^  2 

$(f)  s  S(a,a  )  =  G2(z  ,d)  Ao(z  )  exp[-(-^  + - r—  +  %•)] 

d  0  0  0  o  £  °A^ 

0  0 

2  2 

2a  +a  -2aa  2 

3  S„(z„,d)  exp[ - - ~] 


Equation  (15)  can  be  simplified  by  simple  manipulation  as: 

a.2(l-k  )  (o-k  a .)2 

s(f)  »  S(«)  =  [S  (z  ,d)  exp[-  *  — ]  ]  exp[-  -?  ^-a  ] 

O  kn 

0  0 


(16) 
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where  k  =  — y- - =■  and  a  and  f  are  linearly  related  as  in  Eq.  (4). 

0  «o 

As  seen  from  Eq.  (16),  the  shape  of  the  spectrum  is  Gaussian  but 
shifted  by  an  amount  f  where  f  is 

a  g 


f.  =  -  7  V  a  (-/ - y) 

3  X  3  a  S2o/ 
0  A 


(17) 


This  is  the  single  scatter  value,  as  in  Eq.  (4),  reduced  by  the  factor  in 

brackets.  For  a  given  wind  speed,  the  narrower  A(a)  is,  the  smaller  is 

the  slope  of  f  vs  a,.  The  Gaussian  model  for  A(a),  however,  yields  a 

linear  f  vs.  a,  relationship  that  agrees  well  only  for  small  antenna 
3  3 

pointing  angles  (Birkemeier,  1968). 


Functional  Representation  of  Doppler  Spectrum 
When  Clock-offset  is  Present 

As  mentioned  in  an  earlier  report  that  when  the  transmitting  clock 
(wc)  and  receiver  clock  (wQ)  are  not  identical  then  the  received  Doppler 
spectrum  gets  shifted  by  an  amount  equal  to  clock  offset  Aw  (=wc*w0)-  S° 
if  Act  =  (-4—)  Aw,  then  the  shifted  spectrum  S ( f )  can  be  obtained  from 
Eq.  (16)  as 

(a-Aa-k  a  )2 

S(f)  =  S{a)  =  S0(z0,d,aa)  exp[ - 2  °  3 . ]  (18) 

a  k 
0  0 

2  o2  a  2 

where  S0(z0,d,aa)  =  S0(zQ,d)  exp[-(^-) (- A  y)] 

o  oQ  +2oa 


The  received  spectrum  is  still  Gaussian  but  shifted  by  an  amount  f 
where  f  is: 


f  =  f.  +  2 IT  Aw 

a 


Equation  (18)  can  be  written  in  frequency  domain  as: 


r  (f'ml)  , 

S(f )  =  S0<fa,z0,d)  exp[-  -V—] 


2  2  , 
2  .  1  /°o°A  ,t/2v,2 


2mt 


where  m-|  =  f  ,  m2  =  £(-; ^ - -jlK”) 


V2oA 


(19) 


(20) 


o^+o  2  f  c 

and  Vfa’zo’d)  =  So(zo’d)  exPC"(X7?>(27)  (72,] 

°o  2°A  0 

The  resulting  spectra  is  still  Gaussian.  It  is  well  known  that 
symmetrical  Gaussian  is  represented  by  two  moments  of  the  spectrum  where 
the  first  moment  (m^)  gives  the  centroid  (mean)  and  is  defined  by: 


m 


1 


;fs(i 

f)df 

;s(i 

f)df 

(21) 


and  the  2nd  central  moment  (mg)  gives  the  variance  of  the  spectrum  and 
is  defined  by 

;(f-m1)2S(f)df 
m2  =  /S(f )df 


(22) 


d.u 


The  expression  of  Doppler  spectrum  in  Eq.  (20)  was  derived  based 
on  the  assumption  that  combined  antenna  pattern  and  angular  dependence 
is  Gaussian.  This  results  in  symmetrical  Doppler  spectra  with  mean  m^ 
and  variance  n^.  However,  experimental  evidence  (Birkemeier,  1968)  shows 
that  the  spectrum  tended  to  be  asymmetrical ,  when  the  transmitting  and 
receiving  antennas  are  not  pointed  at  the  great  circle.  At  small  antenna 
offsets  from  boresight,  few  correction  terms  in  Gram-charlier  expansion 
of  the  nearly  Gaussian  spectrum  can  account  for  this  asymmetrici ty  factor, 
e.g.. 


/tTnm, 


i  i  f-m.  2  s  j.  f-m,  .  f-m, 

S(f)  •  gr 


213! 3! 


H6(' 


f-m. 


L)] 


(23) 


Equation  (23)  is  the  corrected  form  of  Eq.  (20)  such  that  total 
input  power  [/S ( f ) df  =  1]  is  unity.  In  Eq.  (23)  is  defined  as  skew 
coefficient  or  asymmetry  coefficient  by: 


where  =  (^) 
excess ,  e.g. 


sk  = 


1.5 


(24) 


_  in- 

(u-rrii )  JS(io)du 

/S(u)du)  and 


is  defined  as  coefficient  of 


Ex  =  4  -  3 


(25) 


where  m. 


,  4  /(w-m. )  s(u)du) 

i  1  \  _ 

'2ti  /  s((jj)du 


For  normal  spectral  slope  Ex=0.  For  other  spectra  having  the  same 
variance  as  the  normal  spectra,  Ex  will  indicate  whether  the  distribution 
has  a  flatter  (Ex^O)  or  sharper  (Ex>0)  peak  than  the  normal  spectrum. 
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One  should  note  that  skew  (sfc)  and  excess  (Ex) can  be  obtained  from 
the  third  (m^)  and  fourth  (m^)  central  moments  of  the  Doppler  spectrum 
If  non-Gaussian  angular  dependence  (£q.  12)  is  considered  and  thus  related 
to  tropolink  parameters.  In  Eq.  (23)  hR(y)  is  a  Hermite  polynomial  of 


degree  n  defined  by: 


H„(y)  =  l-l)n  ey  12  {**  /2) 

"  dyn 


(26) 


with  H3(y)  =  y3-3y  and  H4(y)  =  y4-6y2  +  3  and  Hn+1(y)  =  y  Hn(y)-nHnl (y ) . 


Autocorrelation  Function  of  Doppler  Spectrum 

The  autocorrelation  function  R(t)  of  the  spectrum  density  (Doppler 
spectrum)  S(w)  is  defined  as: 


R(t) 


/  S(u)  e3<i)Tdw 


(27) 


If  S(f)  is  as  given  as  in  Eq.  (23)  then  R(t)  can  be  obtained  by  simple 
manipulation  as 


mV 

R(x)  =  (exp(jm^ t  ))(1- 


3.3 


skm2T 


,c  ,  4  4 
(Ex)m?t 


(28) 


m-j,  m^,  s^.  Ex  are  the  moments  of  the  Doppler  spectrum  S(f)  defined 
earlier  (Eqs.  21,  22,  24,  25). 

The  first  and  second  derivatives  of  Eq.  (28),  which  will  later  be 
used  in  alignment  accuracy  analysis,  is  given  as: 
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2  3  3  r  \  4 

t2  o  9  VEXylTl-  *> 

A(t)  =  (exp(  jm^ t  --^-))  ~T~  T  +^"^_6 —  +  “6~*  T 

(Exfr-m,  skmp  .  (Ex)m~  s.rru  , 

♦J'-JT1*  -¥>- ’■  l-5r*-TT>  ’ 


2  6  2  8 
n1skm2  6  ,  skm2  7- 


72 


T  + 


72 


(29) 


and  R(t)  =  (-B^+B2T^-B2T^+B(jT^-BgT®)-j(c^T+C2T^-C2T^+c^T^) 


(30) 


where 


?  ? 


4  3  (E*)m2 

B2  “  m2  +  mlskm2  +  2 


R  .  9  ,  ,  6  5  2  6  (Ex)mH  ml  skm2 
B3  -  24  (Ex)m2  +  ^  skm2  +  — 24—  +  ■ — 3— 


2  6  8 
„  .  13  2_8  .  mlskm2  .  (Ex)ln2 
“4  72  skm2  72  24 


Br  « 


2  10 
skm2 


5  72 


3  2 

^  2m^ni2 

r  _  m1skm2  .  2skm2  ,  skm2  ,  mlSkm2 

4  t~  +  ~ +  ~r +  — 6 

r  mlskm2  ,  mi(E*)m2  .  skm2 

C3 - 6  +  TT~  +  — 6 


2  8 

r  _  mlSkm2 
^4  36~ 


Accuracy  Analysis  of  the  Slope  Method 


In  the  earlier  report  (Nov.,  1977).  the  error  analysis  of  the  slope 
of  the  quadrature  correlation  function  at  zero  lag  has  been  shown.  The 
formulas  are  summarized  below  for  two  estimators  of  the  slope: 


a)  True  Estimator:  In  this  estimator,  the  slope  (S)  is  estimated  from 
the  expression: 


S  = 


1/  B(t)  £-a(t)dt 


(31) 


with  E [S3  =  Rag(0) 

where  ' E ‘  means  "the  expected  value  of"  and  a(t)  and  e(t)  are  the 
quadrature  components  of  the  complex  process  (a(t)+jB(t))  whose 
spectral  density  is  given  by  Eq.  (23).  The  variance  of  this 
estimator  is  given  by: 


(32) 


where  Raa(x)  =  Re[R(t)]  and  RoB(t)  =  Im[R(x)]. 


b)  Estimator  Hard  Limiting:  In  this  estimator,  the  input  signals  are 

+1  x(t)>0 

first  clipped  by  a  'sqn'  operator  [  0  x(t)=0  ,  where  x(t)  is  the 

-1  x ( t ) <0 

input  signal]  and  then  one  of  the  clipped  signals  is  differentiated, 
e.g.,  the  slope  of  the  quadrature  correlation  function  at  zero  lag 
is  estimated  from: 
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Equations  28,  29  and  30  were  substituted  in  Eqs.  (32)  and  (34)  and  the 
variance  is  compared  for  both  estimators  with  different,  spectral  moments 
(mj ,  m2,  m3  and  m4  is  assumed  zero  for  simplicity).  The  results  are 
summarized  below: 


Effect  of  Fade  Rate  on  Skew:  As  mentioned  before,  the  spectral  moments 
are  functions  of  tropolink  parameters  such  as  wind  speed,  antenna  beam 
width,  angular  dependence,  transmitted  frequency,  the  distance  between  the 
transmitter  and  the  receiver  etc.  The  antenna  elevation  angle  is  assumed  zero. 
The  variations  in  the  link  parameters  causes  these  coupled  spectral  moments 
to  change  accordingly.  So,  simply  varying  these  moments  will  mimic  the 
different  link  conditions,  the  net  effect  of  which  will  result  in  changed 
moments.  In  table  2,  with  m^  =  0,  the  effect  of  fade  rate  under  different 
'skew*  (asymmetric)  conditions  is  shown.  Although  with  different  'skew' 
conditions,  m^  should  correspondingly  be  non-zero,  but  in  calculations 
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first  is  assumed  zero  (Table  2)  and  then  for  a  given  m2,  m^  is  varied 
(Table  3).  For  a  given  skew  coefficient,  the  difference  in  n>]  can  account 
for  the  'residual'  clock-offset  (see  Cq.  19).  From  Table  2,  it  could  be 
observed  that  the  variance  of  both  estimators  increases  in  proportion  to 
the  spectral  width.  If  the  spectral  width  (or  the  fade  rate)  increases 
by  a  factor  of  2,  the  variance  also  increases  by  the  same  factor.  The 


theoretical  explanation  is  given  in  Appendix  1.  The  clipping  loss 

(defined  as  the  ratio  of  ^^^f^^pid^imaT^TSfii^  }  15  ln 
the  range  of  -4  to  -5  db.  If  the  large  amount  of  clock-offset  is  present, 


)  is  in 


then  the  clipping  loss  decreases  accordingly  (see  Table  3).  If  the  Doppler 
spectrum  is  symmetrical  (m^  =  0),  then  it  is  shown  in  Table  4  that  for 
large  clock-offset  (as  increases)  the  clipping  loss  could  result  in 
clipping  gain. 


H 
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Path  UW  Link  RADC  Link 

Antenna  parameters  Symmetrical ,  28’  steerable  Symmetrical,  28'  steerable 

diameter  paraboloids  with  diameter  paraboloids  with 
circular  nominal  -3  db  circular  nominal  -3  db 

bandwidth  of  2.8°,  0°  .  bandwidth  of  .54°  0°  elevation 

elevation  oQ2=2. 82  (degree*1)  oQ2=.105  (degree2) 

Path  length  (2d)  231  km  168  miles,  268.80  km 

Frequency,  wavelength  810  MHz,  .37  m  4.86  GHz , .0625  m 

Size  of  common  volume  110. *5. 6*2. 8  34.9*1 .26*. 63 

(2Lx*2Ly*Lz) 

(lengths  in  km) 

Z  grazing  ray  height  1.57  km  2.12  km 

u  from  link  axis 

Zc,  height  of  earth  .78  km  1.06  km 

surface  from  link-axis 

Zup,  height  of  upper  beam  4.39  km 
from  link  axis 

Zb,  height  of  the  inter-  2.31  km 
section  of  the  trans¬ 
mitter  upper  beam  and 
receiver  from  link  axis 

TQt  single  layer  height  2.5  km  2.5  km 

assumed  from  link  axis 

2  2 

aA  ,  variance  (degree  )  .045, .2 , .547,1 .394 

of  Gaussian  angular  de-  4.905 
pendence  for  fade  rate 
(nu)  of  2, 4, 6, 8  &  10  Hz 
ana  wind  speed  2.5  m/sec 

8,  angle  of  antenna  cen-  1.55°  1.81° 

tral  beam  to  the  link  axis 


Table  D.l 


.001, .005, .014, .032, .077 
(and  .135  for  ni2=  1 1  Hz) 


2.75  km 
2.40  km 


Link  Parameters  for  UW  and  RADC  Links 
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Append!  xQl 


In  this  appendix  the  effect  of  the  change  In  the  width  of  the 
Doppler  spectrum  on  the  variance  of  the  two  estimators  of  the  slope  of 
the  quadrature  correlation  function  at  zero  lag  Is  considered.  Since  the 
variance  is  Inversely  proportional  to  the  averaging  time  and  the  spectral 
width  Is  an  Indicator  of  fade  rate,  so  the  effect  of  the  different  fade 
rates  (due  to  different  wind  velocity,  antenna  patterns,  etc.)  on  the 
required  averaging  time  Is  estimated  by  assuming  that  the  total  received 
power  Is  the  same  under  various  fading  conditions,  e.g.,  if  the  Doppler 
spectrum  Increases  In  width  by  a  factor  of  k,  then 


and  (Al-1) 

R«a<T>  -  R«o<kT> 

where  R^Jt)  Is  the  quadrature  cross-correlation  function  and  Roa(0  Is 
the  quadrature  auto-correlation  function  of  quadrature  processes  a(t) 
and  B(t).  In  deriving  Eq.  (Al-1),  the  well-known  Fourier  transform  scaling 
relationship  Is  utilized,  e.g. 

Substituting  Eq.  (A-l)  In  equations  (32)  and  (34),  one  obtains 


var[S] 


o-  Fr)[R;;(n)"w(,,)*t";«(,’)l2]dn 


and 
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v*r[SCL]  ■  —*•  -  \  [sin  \p(T))]Z 

CL  4TZ  i  (Al-3) 


For  large  ‘T*,  Eq.  IA1-2)  and  (Al-3)  reduce  to 

Var[$]  -  %■/  [R-(n)Rw(n)+[R-e(n)]Z]dn 
and  (Al-4) 

Var[S CL]  -  £  /  dn 

CL  T  0  1-p  -it 

For  Eq.  (Al-4)  It  Is  evident  that  the  variance  of  both  estimators 
Increase  by  the  same  factor  as  the  width  of  the  Doppler  spectrum.  In 
other  words,  an  Increase  In  the  width  of  the  Doppler  rpectrum  requires 
the  Increase  In  the  averaging  time  by  the  same  factor. 
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part  II 


Estimation  of  Quadrature  Correlation  Function  for  all  Delays  by  Analog 
and  Polarity-Coincidence  Method: 


In  the  report  (9/11/77-10/11/77),  the  expected  value  and  the  variance 
of  the  estimator  of  the  quadrature  correlation  function  by  analog  method 
have  been  derived  which  are  summarized  below. 

If  a(t)  and  8(t)  are  the  quadrature  processes  of  unit  power  then  the 
quadrature  correlation  function  Rae(0  is  estimated  by: 

l  T 

Rae(x)  =  y  /  a(t)e(t-T)dt 


(1) 


the  expected  value  and  variance  of  R  .(t)  is  given  by: 

ap 


(2) 


and 

var[Rae(T)j  «  I  }  (l--^)  [R^a(n)+Rae(n+T)Reo(n-T)]dn  (3) 

For  t—0 ,  zero  delay  (lag),  Eq.  (2)  and  (3)  reduce  to: 

ECRae(0)]  =  0  (4) 


and 


where  sgn  is  the  'hard-clipped1  operator,  defined  as 

1  x  >0 

sgn(x)  =  0  x*  0 

-1  x  <0  • 

the  expected  value  of  rc( t)  is  given  by  Van  Vleck  £1966],  e.g. 

E[?c(t)]  a  |  arc  sin(Rae(t)) 

In  the  above  mentioned  report,  the  variance  of  rc(0)  had  been  derived  as 

VartrJO)]  *  -4-  /  [{arc  sinp(n)}2  -  {arc  sinr(n)}2]dn 
c  nzT  -T 

In  Eq.  (8),  p(n)  is  the  autocorrelation  function  of  in-phase  or 
quadrature  processes  and  r(n)  Is  the  normalized  quadrature  correlation 
function  given  by  Eq.  (1). 

In  this  report,  the  variance  of  the  estimator  of  the  polarity- 
coincidence  quadrature  correlation  function  is  computed  for  all  delays 
and  the  two  estimators  of  correlation  function  are  compared  on  the  basis 
of  signal  to  noise  ratio  (SNR).  If  SNR  is  defined  as 

2 

SNR  «  (DC  output  of  the  correlation  at  a  given  lag) 

(variance  of  the  output  at  that  lag) 

then  clipping  loss  or  gain  [Cheng,  1968]  is  defined  as 

(SNR  of  polarity  coincidence  estimator) 

(SNR  of  analog  estimator)  ~ 

Variance  of  Polarity  Coincidence  Quadrature  Correlator 


In  this  correlator,  the  quadrature  correlation  function  (rc(x ))  is 
estimated  by  Eq.  (6).  The  variance  of  this  estimator  is  given  by: 
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Var[r  (x)3  *  X  //  E[sgna(u )sgna(v)sgna(u-T)sgnB(v-r)]dudv-E2[r  (T)]  (11) 

c  V  00  c 

In  Eq.  (11),  one  needs  to  obtain  fourth  product  moment  [McFadden,  1958]: 

E[sgna(u)sgna(v)sgne(u-T)sgnB(v-T)], 

where  a  and  8  are  jointly  gausslan  with  zero  mean  and  equal  variance  (a). 

To  evaluate  this  equation,  we  write: 

y1  *  sgna(u) 
y2  =  sgne(u-T) 
y3  =  sgna(v) 
y4  =  sgnB(v-T) 

and 


P4  *  Pt«(u) 2  0,  B(u-t)  2  0,  a(u) 2  0,  B(v-t) 2  0]  . 

Then  It  Is  known  that,  for  example,  [McFadden,  1966]  and  [Childs,  1967]  that 

p4  *  ts  02) 

The  problem  now  is  to  evaluate  P^.  To  do  this,  we  let 
z1  =  o(u)/a 
z2  *  B(u-t)/o 
Z3  *  «(v)/o 
z4  =  B(v-t)/o 

Then  one  can  obtain  (see  report  9/11/77-10/11/77)  the  covariance  matrix  (C) 

of  the  Gaussian  variables  z^ ,  z2>  z^,  z4  with  zero  mean  and  unit  variance. 

If  n  ■  u-v,  then 


1 

r(t) 

p(n) 

r(n+x) 

r{i) 

1 

r(x-n) 

p(n) 

p(n) 

r(t-n) 

1 

r(T) 

r(n+t) 

p(n) 

r(x) 

} 

P11 

p12 

p13 

p14 

p12 

p22 

p23 

p24 

p13 

p23 

p33 

p34 

p14 

p24 

p34 

p44 

It  Is  well  known  that  P4  for  any  covariance  matrix  can  not  be  obtained 


in  closed  form.  In  the  following  discussion  P4  is  computed  by  Plackett's 
method  (1954)  which  reduces  to  closed  form  expression  (Eq.  8)  for  the  zero 
lag  case. 

We  now  make  use  of  the  result  of  Appendix  1,  Eq.  (12)  and  Van  Vleck’s 
(1966)  formula  ECyiyj3  =  f  arc  sinPjj  for  i,j*l, 2,3,4  and  obtain  the 
expression  for  fourth-product  moment,  i.e. 


E[y1y2y3y4]  -  l  +  4  C2Ii+2I2+i3+i4] 


where  Ij,  I2,  I3  and  I4  are  given  in  the  Appendix  1  (equations  Al-15,  16, 
17,  18). 

Substituting  Eq.  (14)  into  Eq.  (11),  one  obtains 


Var(fc(t))  =  /  (1 -■^-)E[y1y2y3y4]dn]-E2[rc(T)] 

«  1  +  [4  7  (1-M)[2I  +2I2+I3+I4]dn]-E2[r  (x)] 
n<r-T  T  c  3  n 


(14) 


(15) 

(16) 
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Since 


T 

/ 

-T 


(l--^-)l3dn  = 


T 

; 

-T 


l4dn 


Eq.  (16)  can  be  simplified  to 


var[rc(r)]  =  1  +  /  (1  -  M)[2I1+2I3+2I3]dn]  -  E2[rc(t)] 


(17) 


Estimation  of  the  Slope  of  the  Quadrature  Correlation  Function  (R(t))  at 
Zero  Delay,  When  the  Maximum  Value  of  R(t)  is  used  as  an  Estimator. 


In  this  section,  the  'slope'  of  the  quadrature  correlation  function  at 
zero  lag  is  estimated  by  the  maximum  value  of  the  quadrature  correlation 
function,  i.e. 


sl°pe  «  Sm  =  Rm(ro)-R(0) 

To 

Here  Rm(t0)  is  the  estimated  maximum  value  of  R(t)  occurred  at  lag  t  . 

If  R(t)  is  estimated  by  the  Analog  method,  (see  Fig.  1),  then  the 
slope  estimator  is  given  by: 


AN 


lyy-wi 


the  expected  value  of  this  estimator  is  given  by 

.W 

T  T 

0  0 

If  R(t)  is  estimated  by  the  Polarity-coincidence  method »  then 


e£san^  = 


;  .  rc<ToK(0) 

spcc  '  r0 

the  expected  value  of  this  estimator  is  given  by 
»  EtrfvJ]  ,  n  i 

— v —  !rrc(To)s!rarcsinyT0) 


(18) 


(19) 


(20) 


(21) 


e£spcc^  = 


(22) 
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Variance  Calculations: 

The  variance  of  the  slope  estimator  is  given  by: 

T*var[Sm]  =  var[Rm(x0)]  +  var[R(0)]  -  2E[Rm(xQ)R(0)]  (23) 

The  variance  of  Rm(xQ)  an<^  R(0)  can  be  computed  by  the  methods 
discussed  in  the  earlier  section  for  both  analog  and  polarity-coincidence 
cases.  Now  a  method  is  derived  to  compute  the  third  term  for  both  cases. 

When  R|11(t0)  and  R(0)  are  computed  by  Analog  method,  then 

1  T 

R(0)  =  f  /  a(t)  e(t)dt 
1  0 
and 

VTo>  *  W  ■  T  '  <24> 

So  E[Rm(To)R(0^  is  given  by: 

!  T  T 

=  "V  Ef'f  «(u)e(u)du  I  a{v)e(v-T  )dv] 

mo  y*  g  0  0 

,  TT 

=  -4  //  E(a(u)e(u)o(v)e(v-x  ))du  dv  (25) 

r  oo 

If  a  and  6  are  jointly  Gaussian,  then 

E[a(u)e(u)a(v)e(v-x0)]  =  Rag(u-v+x0)  RBa(u-v) 

+  Raa(u_v)  Ree(u'v+To)  (26) 

If  u-v  *  n.  then  Eq.  (25)  can  be  rewritten  as 

■  Ety. „)#«»]  =  }  7  (l  -  («aa(n)R68(n«0)  -  V>lV"%))d'1  (27) 

A  A 

If  Rm(r0)  and  R(0)  are  computed  by  polarity  coincidence  method,  then 
1  T 

R(0)  *  r.(0)  =  y  /  Sgn  a  (t)  Sgn  e  (t)  dt 
c  1  0 


(28) 


53 


.32 


and 


R«('o>  * 


rc<V 


/  Sgn  a ( t)  Sgn  e(t-Tft)  dt 
0  0 


(28) 


So  E[rc(To)rc(0)]is  9<ven  by: 


TT 


EtrJiJrJO)]  =  -U //  E[Sgn  a(u)  Sgn  e(u)  Sgn  a(v)  Sgne(v-t  )]  du  dv  (29) 

c  0  c  t  00  0 


If  yy]  =  Sgn  a(u) 
yy2  »  Sgn  B(u) 
yy3  =  Sgn  a(v) 
yy4  =  Sgn  b(v-t0) 

then,  using  Eq.  (12)  and  Eq.  (A3-12),  one  can  obtain  the  expression  for  the 

fourth-product  momment  E[yy^  yy2  yy3  yy43  (similar  to  Eq.  (14)), 

ECyy,  yy2  yy3  yy4]  ■  1  +  Cvi?  ♦  vi2  +  vi3  +  vi4]  (30) 

where  VI,,  VI2>  VI3>  and  VI4  are  given  in  Appendix  3  (Equations  A3-8,9,10,ll). 

Substituting  Eq.  (30)  in  Eq.  (29),  one  obtains 

E[rc(xo)rc(0)]  =  j  /  Jl  --ly^ECyy,  yyg  yy3  yy4]  dn 

»  1  +(^)  }/  (l  -Mj  (VI,  +  VI2  +  VI3  +  VI4)  dn  (31) 


Substitution  of  Eqs.  (27)  and  (31)  in  Eq.  (23),  one  obtains  the  variance 
of  the  slope  estimator  [S  ]  for  both  anlaog  and  polarity-coincidence  methods 
and  using  Eqs.  (3,5,8,17,20,22,23)  and  Eqs.  (9)  and  (10),  one  can  compare  the 
two  estimators  on  the  basis  of  signal  to  noise  ratio. 


t- 
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Functional  Form  of  Doppler  Spectrum 

In  the  report  (12/11/77  -  1/11/78),  the  asymmetrical  Doppler  spectrum, 
at  small  antenna  offsets  from  boresight,  was  represented  by  Gram-Charlier 
expansion.  This  expansion  has  been  modified  here  in  such  a  way  that  the 
first  four  central  moments  of  the  expansion  are  Identical  to  that  of  the 
spectrum  derived  from  the  combined  antenna  pattern  and  angular  dependence. 
The  modified  Gram-Charlier  expansion  of  the  spectrum  S(f)  is  written  as 


3m^ 


The  autocorrelation  function  R(T)  of  the  spectrum  S(f)  defined  as 

R(t)  =  /  S(w)  e^w*d«  (33) 

-•0 

is  obtained  by  manipulation  as: 


(34) 


The  first  derivative  of  Eq.  (34),  which  will  be  used  In  alignment  accuracy 
analysis  Is  given  as: 


nuT 

ft(T)  =  (exp(j  mj  t  -  -7—)) 


<  J*5 

.  m  3  3v  2 

J -j(  ^  "  2  )T 


+  ( 


15  2  A 

ml s Km2 *  ,  ^Ex^m2  mV 


)t” 


,c  .  2  2.5  <r  .  3  2  3 

/(Ex)m2m1  _  s^  %  4  /(Ex)m2  ^  s^  5 

24  6  'T  *  '  24  12  ;x 


2  3  2  4 

in^S^  6  sKm2  7 

j  72  T  +  -72“t 


(35) 


Accuracy  Analysis  of  the  Correlation  Method  of  Alignment 

In  the  report  12/11/77-1/11/78,  it  was  shown  that  the  spectral  moments 
are  functions  of  tropol ink-parameters  such  as  wind-speed,  antenna  beam 
width,  angular  dependence,  transmitted  frequency,  the  distance  between 
the  transmitter  and  the  receiver,  etc.  The  antenna  elevation  angle  is 
assumed  zero.  The  variation  in  the  link  parameters  cause  these  coupled 
spectral  moments  to  change  accordingly.  So  simply  varying  these  moments 

1 

will  mimic  the  different  link  conditions,  the  net  effect  of  which  will  result 
in  changed  moments.  In  Table  1,  with  m^=0,  the  effect  of  fade  rate  under 
different  'skew'  (asymnetric)  conditions  Is  shown.  Although  with  different 
'skew'  conditions,  m^  should  correspondingly  be  non-zero,  but  In  calculations, 
first  m1  is  assumed  zero  (Table  1)  end  then  for  a  given  n^,  is  varied 
(Table  2).  For  a  given  skew  coefficient,  the  difference  in  mi  can  account 
for  the  'residual  clock-offset'.  From  Table  1,  It  could  be  observed  that 
if  the  spectral  width  (or  the  fade  rate)  Increases  by  a  factor  of  2,  the 
variance  decreases  by  the  same  factor.  The  theoretical  explanation  is  given 


In  Appendix  2.  The  clipping  loss,  defined  as  In  Eq.  (10),  Is  In  the  range 
of  -1.5  to  -3  db.  In  Table  3,  the  case  of  the  symnetrlcal  Doppler  spectrum 
Is  considered.  In  Table  2  and  3,  It  should  be  noted  that  In  some  cases 
the  clipping  loss  could  result  In  the  clipping  gain.  It  would  be  of 
significant  Interest  to  Identify  those  cases  a  priori.  From  Tablesl,  2,  and 
3  It  can  be  concluded  that  the  variance  decreases  by  the  same  factor  as 
the  averaging  time  Is  Increased. 

If  the  slope  Is  estimated  from  the  correlation  function  (Eq.  18),  then 
it  could  be  observed  from  the  equations  23,  27  and  31,  that  the  variance  of 
the  slope  In  both  cases  will  decrease  as  the  spectral  width  (or  fade  rate) 

Is  Increased  which  Is  different  than  that  of  the  slope  estimation  by  True 
and  hard-limiter  method  (report  11/12/77-1/11/78)  where  the  variance 
Increases  In  proportion  to  the  fade  rate.  So  It  can  be  concluded  that 
the  slope  estimation  by  Eq.  (18)  Is  'better'  than  the  estimation  by  True 
and  hard-limiter  methods. 


D.  36 


References 

1.  M.  C.  Cheng,  "The  clipping  loss  In  correlation  detectors  for 
arbitrary  Input  signal  to  noise  ratio,"  IEEE  Trans,  on  Infonaatlon 
Theory.  Vol.  IT- 14.  No.  3.  Nay  1968,  pp.  382-389. 

2.  D.  R.  Childs,  "Reduction  of  the  nultl variate  normal  Integral  to 
characteristic  form,"  Blometrlka  54,  293-300,  1967. 

3.  J.  A.  Mcfadden,  "The  fourth  product  moment  of  Infinitely  clipped  noise," 
IEEE  Trans,  on  Information  Theory  4,  159-162,  1958. 

4.  R.  L.  Plackett,  "A  reduction  formula  for  normal  multivariate  Integrals," 
Blometrlka,  Vol.  41.  pp.  351-369,  1954. 

5.  J.  H.  Van  Vleck  and  0.  Middleton,  "The  spectrum  of  clipped  noise,"  Proc. 
IEEE,  Vol.  54.  pp.  2-19,  1966. 


D.37 


D.  38 


! 


where 


r  -  r  -  1-a2-d2-b2+2abd 

11  44  “jPT 

.  (a-a3-bd*acd+ab2-bc) 

12  34  jPl 


c13  *  c24 


ad-a2b-b+ac+b3-bcd 


.  a2d-2ab+c+b2d-cd2 

14  '  I?1! 


C22  *  C33 


l-a2*b2-c2+2abc 

ic^l 


C23  "  ' 


d-2ab+a2c+b2c-c2d 

I?1] 


and  |C_1|  «  1-2a2-d2+4abd-2b2+a4-2a2cd-2a2b2+4abc+b4-2b2cd-c2+c2d2. 

The  method  followed  to  compute  P4  is  the  application  of  the  general 
method  given  by  Plackett  [1954],  Let  the  correlation  coefficients  be 
denoted  by  p^j  (Eq.  Al-2),  then  since  the  density  function  obeys  the 
relation: 


a* 

ap 


we  have,  for  example 

aP, 


ap 


12 


U 

00 


mm  2 

//*-* 
00 


a2a 


1j  3xiaxj 


3x^X2 


dx^dx  ■]  dx3dx4 


1  P2<c33*c34,c44) 

2n(l-p12  ) 


(Al-5) 


(Al-6) 


There  will  be  6  such  equations,  one  for  each  correlation  coefficient. 
Equation  (Al-6)  Is  a  special  case  of  Eq.  (6)  given  by  Plackett  [1954]. 
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Since 


P2tc33*c34’c44l  “  }  +  21  arc  s1n 
combining  Eq.  (Al-6)  and  (Al-7),  we  have 


-c 


34 


9P. 


1 


9p12  8n(l-p^2) 


1  *  -  arc  sin 

n 


(c33c44)a 


L  (c33c44] 


There  are  six  such  equations,  upon  combining  them  we  obtain 


p«  ■  W  *  b 1  / 


1J 


ij  L 


1  +  |  arc  sin 


(cppcqq)2 


JLNj. 


where  (p.q.i.j)  Is  a  permutation  of  (1,2, 3, 4). 

If  P^kg]  Is  the  normal  orthant  probability  of  four  Gaussian  variables 
of  zero  mean,  unit  variance  and  covariance  matrix 


’ 1  0  1  O' 

1  0  1 
1  0 
1 


then 


at  •  *  -  ^ 


(Al-7) 


(Al-8) 


(Al-9) 


( AT  -  TO) 


m 


'«  ■  VM  *  n  *iirc  sin  -(-V  J  - : 12 


^C33C44)  /\-)£ 


12 


^  O  ("Ca*)  d^l a  . 

+  2  /  [1  + 1  arc  sin  - r]  - — 


^c22c44^  A-X* 


13 


2  ^"C23^ 


+  /  [1  +  £-  arc  sin 


+  /  [1 +  f  arc  sin 


dx 


14 


(c22c33)*  /_x2 


(-c14)  dx 


23 


(ciic44)  A-x\ 


14 


■3 


23 


(Al-11) 


P4[kg]  can  be  computed  from  Cheng  [1968]  e.g. 

p4[k03  =  of  Cheng  =  ^ 

Combining  (Al-11)  and  (Al-12),  one  can  obtain 


(Al-12) 


P4  =  F  +  inf  C2sin"1a  +  2sin"^b  +  sin_1c  +  sin_1d]  +  — 


FPM 


where 


(Al-13) 


a  -c 

FPM=[  2  /  arc  sin[ 

0 


”^34  dx.p  b  dx., 

—  -  ■  J  -^J+V{*rc  sin[—  ]  -JJj 

(c33c44>  /l-x^2  (C22c44*  /Tx^ 

+  [/  arc  sin[ - — r]  — =•■*  3+f7  arc  sin  — --  ^t] 

0  ic22c33^  /Tx^  0  *cllc44^  \Znx23 


(Al-14) 
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Substitution  of  from  (Al-5)  into  (Al-14),  one  obtains: 


FPM  =  21,  +  212  +  I3  +  I4 


where 


I  .  7-48- 


arc  sin 


_ (A-A3-BD+ACD+AB2-BC) _ 

( 1 -A2-B2-D2+2ABD  f (1 -A2-B2-C2+2ABC f 


i 

2  1  n 


'l-B 


arc  sin 


(B+A2B-AD-AC-B3+BCD) _ 

( 1  - A2 -  B2  -  D2 + 2  ABD  ( 1  -A2-B2-C2+2ABC)* 


l3Sf 


dC 


and 


0 


arc  sin 


l,  *  / 


dD 


4 


arc  sin 


(D-2AB+A2C+B2C-C2D) 
( 1 - A2- B2- C2+2ABC ) 


(A20-2AB+C+B2D-CD2) 
( 1 - A2-B2-02+2ABD ) 


Since  A,  B,  C,  0  are  interrelated  variables  as  in  Eq.  (Al-3),  so  variation  in 
one  results  in  the  variation  of  another.  So  the  above  integrals  I,,  1^. 

I3,  I ^  can  be  changed  as: 


I  «  /  r(t’  )dif_.  arc  s,n 
0 


I.  ■  ",  .re  sin 


(r(t  * )-  r?(t  *  )-br(t  ’-n)+r(t  ’  )r(t  ’+n)r(t  ’-n) 
+r[t_‘ J+ b  -br(t  *+n) ) 


^7(n‘) 


(l-r2(t,)-b2-r2(t’-n)+2r(t,)(b)r(t,-n))V 
(l-r2(t,)-b2-r2(t*+n)+2r(t,)(b)r(t,+r))% 
(p(n')+a2p(n,)-a(r(T-n’)+r(T+n,))-p3(n,r| 

+p(n')r(T+n')r(T-n') _ 

(l-a^-p2(n,)-r2(x-n')+2ap(n,)r(T-n,))% 
{l-a‘’-P2(n,)-r2(T+n')+2ap(n')r{T+n'))*  J 


(Al-15) 


(Al-16) 
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Appendix  2 

In  this  Appendix,  the  effect  of  the  change  in  the  width  of  the  Doppler 
spectrum  on  the  variance  of  the  two  estimators  of  the  quadrature  correlation 
function  at  all  lags  is  considered.  Since  the  variance  is  inversely 
proportional  to  the  averaging  time  and  the  spectral  width  is  an  indicator  of 
fade  rate,  so  the  effect  of  the  different  fade  rates  (due  to  different  wind 
velocities,  antenna  patterns,  etc.)  on  the  required  averaging  time  is  estimated 
by  assuming  that  the  total  received  power  is  the  same  under  various  fading 
conditions,  e.g.  if  the  Doppler  spectrum  increases  in  width  by  a  factor  of 
k,  then 


-  RaB(kT) 

and  (A2-1 ) 

R  (x)  -  R  (kx) 

act  aa 

where  Rag(x)  is  the  quadrature  cross-correlation  function  and  R^-r)  is  the 
quadrature  auto-correlation  function  of  quadrature  processes  a(t)  and  B(t). 

In  deriving  Eq.  (A2-1),  the  well  known  Fourier  transform  scaling  relationship 
is  utilized,  e.g. 

iHjl-flax) 

Substituting  Eq.  (A2-1)  into  Eqs.  (3)  and  (15),  one  obtains,  for  large  T, 

Var[Rae(x)]  =  ^  ^  (R^(n)+Ro6(n+T)Rga(n-T)dn  (A2-2) 

and 

,  T 

Var[rc  (t)]  *  jj  /  Efy^^y^dn  (A2-3) 

where  rc(t)  is  the  clipped  quadrature  correlation  function. 


From  Eqs.  (A2-2  and  A2-3),  it  is  evident  that  the  variance  of  both 
estimators  decreases  by  the  same  factor  as  the  width  of  the  Doppler  spectrum 
increases.  In  other  words,  an  increase  in  the  width  of  the  Doppler  spectrum 
requires  the  decrease  in  the  averaging  time  by  the  same  factor. 
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Appendix  3 

In  this  Appendix,  the  expression  for  (quadri variate  normal  orthant 
probability)  is  derived  which  is  required  in  computing  the  variance  of  the 
slope  estimator  by  polarity  coincidence  method.  The  method  is  described 
in  detail  in  Appendix  1  and  has  been  applied  in  comouting  the  variance  of  the 
quadration  correlation  function  estimator  by  polarity  coincidence  method. 

In  the  case  under  discussion,  the  covariance  matrix  (c)  is  given  by 


'1  0  a  b  ' 

1  c  d 

c  = 

l  g 

l 

where  a  =  p  (n ) 

b  =  r(n+TQ) 

c  «  -r(n)  *  r(-n) 

d  =  p(n+TQ) 

g  =  r(t0) 

The  inverse  of  the  covariance  matrix  (c)  is  given  by: 


C11  0  c13  c14 

-1  ^  c22  c23  c24 

c33  c34 

c44 


(A3-1 ) 


(A3-2) 


(A3-3) 


'11 


_  l-q^-c^+2cdq-d2 


|c 


-1 


where 
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-a+hq+ad  -bed 

Ci<j  “  “  31 

13  I?1! 

2 

ag-b-acd+bc 
C’*‘  |c-'r| 

.  l-a2-D?-g2*2at>q 
2 

c-gd+abd-b  c 

C23=‘  k'1! 

2 

cq-d+a  d-abc 

‘24  =  1C1! 

i-b2-d2 
33  =  Ic-1! 


'34 


q-cd-ab 

Ic"1! 

1-a2-c2 


44  i-l 


Ic"1  I 


and 


jc_1|  =  i-g2-c2+2cgd-d2+(-a2+2abg+a2d2-2abcd)-b2+b2c 


,2.u2„2 


Similar  to  Eq.  (Al-13),  one  can  write 


p.  =  1  +  J-  [arc  sina+arc  sinb+arc  sinc+arc  sind+arc  sing]  +  — v 
4  8  8n  L  4n 


FPM 


where 

FPM 


a  -c,.  dx,,  b  .  "c23  ,  ^14  i 

~  ^  *  [,o m  s,"‘ W1  ^ 

'  ,_a13  14 


+  [/  arc  sin[ 
0 


inCT^r]  ^  +  E(  arC  SinC{F^3 


(A3-4) 


(A3-5) 


(A3-6) 
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Substituting  from  (A3-4)  into  (A3-6),  one  obtains 


FPM  =  VIj  +  V I2  +  VI3  +  VI4 


(A3-7) 


where 


VI,  =  / 


dA 


arc  sin 


1  /-A2 


(D+ABC-CG-A^D) _ 

( 1 - A2-B2 -G2+2ABG  )*  ( 1  - A2- CZ  ? 


vi2  =  / 


dB 


arc  sin 


(C-GD+ABO-BC) _ 

( ) -A2-B2-G2+2ABG)% ( 1 -B2-D2)4 


VI,  =  / 


dC 


0  A-CZ 


arc  sin 


(B+ACD-AG-B(r ) 


(1 -G2-C2-D2+2CDG)% ( 1 -A2-C2 )% 


VI4  =  1 


dO 


1  J\- O2" 


arc  sin 


(A+BCD-BG-AD^) 


{ 1 -G2-C2-D2+2CDG)* (1  -B2-D2)% 


Since  A,B,C,D  are  interrelated  variables  as  in  Eq.  (A3-2),  the  variation 
in  one  results  in  the  variation  of  another.  So  the  above  integrals  Ij, 

I2,I3,  and  I4  can  be  changed  as 

(p( n'+Tft)-p(n’)r(n’+T  )r(n,)+qr(n')-p2(n,)p(n,+T0) 


VI,  s  f  —  arc  sin 

1 


( (1-P2(n' )-r2(n'+Tn)-g2+2gp{n' )r(n‘+x  ))4 


(l-p2(n,)-r2(n,))%  ) 


(A3-8) 


VI,  =  /  arc  sin 


0  /T72(f) 


( c-r(t,-n)p(t,)+ar(t,)p(t,)-rz(t'Jcl. 


|((  1  -a^- r2 ( t ‘  )-r2(t'-n)+2ar(t’  )r(t’  -n))^ 
(l-r2(t')-p2(t'))^) 


.+  °j  ^(nl  )-dn —  arc  sin 


0  /T-r^Cn* ) 


(-rfn'-x.l-gpln' )+p(n'-T  )r(n' )p(n‘ )+r2(n' )r(n'-t  J 


j((l-pi:(n,-T0)-r'i:(n,)-g?+2p(n,-T0)r(n,)g) 
(l-r2(n‘ )-P2(n' ))%) 


(A3-9) 


♦n  r(n<\(in'  («*(n *+t )-p(n ' )r(n* )p(n'+t  )-p(n ' )g-rc(n 1  )r(n '+t J 

VI  s  .  /  rin  JQn  arc  sin  _ 2 _ * _ 2 _ 2_ 

3  0  <CV"j  ((lV-r^n^-p^n'+T^-egrtn'Jpfn1^))* 


(1-p2{n*)-r2(n‘))*) 


(A3- 10) 


VI-  =  °f  ■  arc  sin 

0  /-p2(f) 


(a»r(tl)cp(t')-r(tl)r(t‘-n)-ap2(t1)) 
((l-c2-P2(t‘ )-r2(t'-n)+2cp(t‘ )r(t'-n) )* 
(l-r2(f)-p2(f))*) 


Vn  P(n* )<jn‘  (<,^n,"T0^_r(n' )r(n,*T0)p(n' )-r(n' )g-p(n'«To)p2(n' )) 

+  /  — * — L - arc  sin  - s — g - « - ; - 


0  ^-P2(n') 


({l-g2-r2(n,-To)-p2(n')-2gr(n,-To)p(n')  )* 
(l-r2(n*)-p2(n'))%) 


(A3-11) 


Substituting  equations  ( A3-8 ,9,10,11 )  iff  Eq.  (A3-5),  one  obtains 


1  2 

P4  =  [1+1+  — [arc  sina  +  sin  arcb  +  arc  sine  +  arc  sind  +  arc  sing] 


+  7  tWW] 


(A3-12) 


At  the  present  there  is  an  increasing  interest  in  the  use 


of  communication  channel  simulators  for  the  checkout  and  comparative 
testing  of  modems.  The  cost  and  difficulty  of  obtaining  operating  links 
and  the  impossibility  of  comparative  modem  testing  with  the  same  link 
(at  the  same  time)  have  been  the  factors  stimulating  this  Interest. 
Another,  less  obvious  factor,  Is  the  need  for  a  simulator  as  a  "standard 
channel"  to  be  used  In  the  calibration  of  channel  measurement  equipment. 
The  only  meaningful  checkout  of  channel  measurement  equipment  (and  also 
modems)  prior  to  on-the-air  use  involves  a  synthetic  channel  which  can 
be  set  up  to  simulate  known  channel  parameters. 

The  mathematical  model  around  which  the  simulators,  In  the  past, 
have  been  built  is  that  of  a  band-limited  uniform  tapped  delay  line 
[Baghdady,  1969.  Bello,  1969 ;  Oaly,  1964.  Fitting,  1967.  Klein,  1967; 
Hancock',  1966;  Freudberg,  1965;  Rose,  1966;  Stein,  1965;  Walker,  1966] 
with  each  tap  undergoing  Independent  fading  according  to  complex 
Gaussian  statistics.  The  choice  of  complex  Gaussian  statistics  can  be 
justified  on  the  grounds  of  central  limit  theorem  as  applied  to  the 
contribution  at  each  path  delay  caused  by  scattering  from  a  large 
number  of  atmospheric  refractive  Index  Inhomogenietles.  Actually,  It 
has  been  noted  [Bennet,  1948]  that  as  few  as  six  sine  waves  with 
Independently  fluctuating  random  phase  will  give  a  fluctuating  result 
whose  envelope  closely  follows  Rayleigh  statistics  and  whose  phase  Is 
uniformly  distributed.  Statistically,  the  problem  of  determining  the 
resultant  of  a  set  of  random  phasors  Is  familiar  as  that  of  two-dimensional 
random  walk.  A  more  basic  view  of  the  two-dimensional  random  walk  follows 
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by  regarding  the  constituent  Independently  varying  phasors  to  be  resolved 
into  quadrature  components.  The  simultaneous  values  of  the  quadrature 
components  of  the  resultant  at  any  instant  are  then  readily  shown  to  be 
uncorrelated,  and,  by  the  central  limit  theorem,  to  Independently  approach 
Gaussian  variates  as  the  number  of  steps  in  the  random  walk  (the  number 
of  contributing  phasors)  becomes  large. 

With  regard  to  experiments,  the  short-term  fading  statistics  of  a 
received  signal  have  often  been  observed  to  have  a  nearly  Rayleigh 
distribution  which  is  consistent  with  the  complex  Gaussian  assumptions. 

It  Is  well  known  that  the  received  signal  is  the  linear  combination  of  the 
signals  from  the  different  paths  (delays).  The  dlstinguishability  of 
several  paths  implies  that  they  are  resolvable,  i.e.,  their  electric 
lengths  are  sufficiently  different.  More  concretely,  the  situation  usually 
described  as  'multipath'  is  one  in  which  more  than  one  propagation  path 
exists,  with  the  time  differentials  between  the  several  paths  larger  than 
the  period  of  the  highest  'significant'  frequency  components  In  the 
modulation  of  the  transmitted  signal. 

Shortcommings  of  the  Past  Simulators.  The  above  mentioned  simulators, 
though,  give  adequate  performance  In  the  laboratory,  but  they  have  several 
shortcommings  from  a  meteorological  and  communication  point  of  view. 

1.  Aligned  antennas  are  assumed.  Normally  antennas  are  less  than 
optimally  aligned,  so  the  Doppler  spectrum  from  the  single  tap 
is  not  necessarily  synmetrlcal  and  the  quadrature  components 
therefore  become  statistically  dependent  processes. 

2.  Zero  clock-offset  Is  assumed.  Normally  there  is  clock-offset 
between  the  transmitting  and  receiving  clocks,  e.g.,  the  clock 
frequencies  are  not  identical.  Also  there  is  some  drift  present 
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In  the  clock,  though  the  clocks  used  are  usually  fairly  stable 
(1  in  1010Hz).  Both  the  effects,  clock-offset  and  drifts  of  the 
clocks,  should  be  included  in  the  effective  channel  simulator. 

3.  There  is  no  provision  to  include  the  average  horizontal  drift  of 
the  scatterers--a  factor  which  is  vitally  Important  to  a  realistic 
model . 

4.  There  is  no  provision  to  include  the  angular  dependence  of  the 
scatterer  which  assigns  a  signal  strength  to  the  scatter  as  a 
function  of  the  scattering  angle. 

5.  There  is  no  provision  to  directly  include  meteorological  conditions. 
These  simulators  are  not  responsive  to  a  change  from  turbulent 
conditions  (fast  fading)  characterized  by  a  Rayleigh  envelope  to 

an  anisotropic  medium  (slow-fading)  usually  characterized  by  a 
Rician  distributed  envelope  except  by  "brute  force"  method,  1,e.t 
forcing  one  scatterer  to  be  very  strong. 

6.  Airplane  free  common  volume  is  assumed.  Normally,  while  link 
is  In  operation,  there  are  airplanes  flying  in  the  common-volume. 
Reflections  from  airplanes  can  be  viewed  as  the  reflections  from 
a  "large  scatter"  resulting  In  the  'v' shape  Instantaneous  frequency 
of  the  return  complex  signal. 

In  the  following  description,  a  channel  model  Is  proposed  in  which 
the  above-mentioned  effects  are  included.  This  discussion  will  deal  with 
the  simulation  of  single  tap  In  the  tapped  delay  line  model  (or  single 
path)  of  the  channel.  This  model,  in  general,  can  be  used  for  other 
channels  (UHF,  VHF,  Ionospheric  channel,  underwater  channel,  etc.).  This 
channel  model  Is  divided  Into  three  parts  as  shown  In  Fig.  1. 
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Models  for  airplane  reflections  will  be  the  topic  of  further  research 
and  will  be  discussed  at  a  later  date.  In  the  following,  a  model  for 
diffuse  scattering  such  as  would  occur  from  a  single  turbulent 
atmospheric  layer  is  discussed,  and  then  a  model  for  low  frequency 
specular  component  is  described. 

Scatter  Model.  As  mentioned  earlier,  the  models  in  the  past  utilized 
the  symmetrical  Doppler  spectrum  with  no  clock-offset,  the  tap-gains 
were  varied  to  represent  the  different  fading  conditions.  In  the  proposed 
model,  the  other  parameters  of  the  tropo-channel  are  also  included  which 
mimic  the  effects  of  the  clock-offset,  misalignment  of  the  antennas, 
spatial  non-statlonarities  In  the  drifting  scattering  layers,  etc.  In 
the  following  discussion,  only  single  tap-delay  line  (single  path)  model 
Is  considered  for  simplistic  reasons. 

In  the  report  12/11/77  -  1/11/78,  It  was  shown  that  the  received 
Doppler  spectrum  S(f)  could  be  represented  as: 


S(f)  -  C0F[GT(f)  A(f)  GR(f)] 


(1) 
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where  G^(f)  and  GR( f )  are  assumed  Gaussian  transmitter  and  receiver 
gain  functions,  A(f)  is  the  Angular  Dependence  function  and  COF  is 
receiver 'clock-offset  shift'  operator.  If  there  is  no  offset  then  received 
spectrum  SQ(f)  can  be  represented  as: 

So(f)  ■  GT(f)  A(f)  GR(f)  (2) 

S(f)  and  5 Q(f)  are  related  by:  SQ{ f)  ■  S(f+af)  where  af  is  the  clock-offset. 

If  the  transmitter  and  reciever  antenna  beam-widths  are  identical  with 
Gaussian  beam  width  "measure"  of  and  the  transmitter  is  assumed  pointed 
at  the  horizon  over  the  great  circle  but  the  receiver  antenna  is  azimuthally 


pointed  at  an  angle  aa,  then  one  obtains: 

GT(f)  “  60  ««P(-  -4 ) 

°o 

(3) 

and  GR(f)  *  G0  exP(~^  ) 

°o 

(4) 

In  equations  (3)  and  (4),  (Merril,  1968) 

■  * » ’  *"d  *5  - 

(5) 

Where  v  is  the  cross-path  wind  velocity,  x  is 

transmitted  CW  signal,  2d  is  the  length  of  the 

the  wave  length  of  the 

link  axis,  a ^  is  -3db 

beam-width  of  a  single  antenna  pattern  and  Gq  is  constant. 

Similarly  angular  dependence  function  is  written  as:  (Merril,  1968) 

A(f)  -  A0[f2*kJ]’11/6  <6> 

Where  AQ  Is  constant  depending  upon  the,cross-path  wind-velocity  and  intrinsic 
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strength  of  the  scatterers  and  kQ  is  given  by: 

„  z°  •  2v 

Ko  "  T  ~  (7) 

20  Is  the  height  of  the  layer  In  consideration.  For  small  antenna 
misalignment  and  narrow  beams,  A(f)  can  be  approximated  by: 

A(f)  -  A0[f2+K^r2  (8) 

Simulation  of  the  Complex  Channel  :  In  this  discussion,  the  tropo-channel 
will  be  simulated  and  then  the  clock-offset  effect  will  be  included  In 
the  overall  simulation.  It  is  well  known  that  a  process  with  spectral 
density  SQ(f)  can  be  generated  by  passing  white  noise  through  the  filter 
of  transfer  function  H(S)  such  that: 

|H(jw)|2  -  S0(f)-ST( f)  A(f)  6R(f)  (9) 

Similarly,  If 

|HT(ju))|2  -  GT(f) 

|HA<  ju.)  I2  -  A(f)  (10) 

and  |HR(  jw)  |2  ■  GR{ f) 

then  one  can  write 

|HT(Jw)|2  -  |HT(j«)|2  |HA(j«)|2  |HR( ju) |2  (11) 

Using  equations  (3),  (4),  (6),  (10)  and  (11),  one  obtains 

2  f2 
|HT(ju)r  -  G0  exp  - 

°o. 


(12) 
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and  •  |H(ju 

eqs.  (12)  to 

White  Noise 

White  Noise  ~ 

White  Noise  ' 

White  Noise  — 

White  Noise  — 


lHR(jw)|  -  %  exp 


Lf-u)* 


|HA(ju)|2  -  Ao  [f2+K2]’2 


'I2  *  Go  Ao  «P 


e) 

°0 


(15)  are  shown  In  block-diagram  form  In  Figure  2. 


(13) 

(14) 

(15) 


SQ(f)  -  |H(Jw)|2  -  GT(f)  A(f)GR(f) 


'clock  offset 
shift'  operator 


S(f) 


ht(J«)  h>6T(f) 


^Jw) 


A  H.(Jw)  f->G.(f) 


^  Hr(J«)  \->G0(f) 


?  f2 

|HT(J«)r  “  Goexp  ( *  -S ) 

°0 

|HA(jw)|2-  Aq  (f 2+K2]'2 

|Hr(»|2  »  Goexp(-  ) 

0 


•* ■  - 1 


i-5 

i 

HT( ju>)  - 

,  H|( Ju) 

—  3 

hr(j«) 

f 

C  _ 

H(j«) 

■  - 

-  -  -  -  • 

s  (f) 

>° 


clock- 

offset  shift 
operator 


S(f) 


Block  diagram  forms  of  equations  12-15 
Fig.  2 
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Some  Properties  of  |H(jg>)|  .  Here  It  Is  shown  that  the  asymmetric 
filter  H(S)  causes  the  asymmetrical  Doppler  spectrum  which  results  In  the 
cross-correlation  function  of  the  output  (of  the  filter)  quadrature 
processes  (inphase  and  quadrature).  It  has  been  mentioned  in  earlier 
reports  that  the  misalignment  of  the  antennas  and  the  clock-offset  both 
result  in  the  output  quadrature  processes  being  correlated.  It  will  be  shown 

/  i 

in  the  later  section  that  by  decomposing  H(S)  in  even  and  odd  components, 
one  can  directly  relate  the  odd  part  of  H(S)  to  the  cross-correlation  of 
the  output  processes. 

2 

The  even  and  odd  parts  of  |H(jw)|  are  defined  as: 

Ev|H(ju)|2  =  Even  part  of  |H(ju)|2  «  (1/2  )|H(  ju>)  |2  +  (1/2  )|H(-joi)  | 2 

(16) 

and  0D|H(ju>)|2  *  Odd  part  of  |H(  ju) |2  «  0/2  )  H(  ju)  |2  -  (1/2  )|H(-ju)  |2 

2 

From  equation  (16),  it  is  clear  that  the  odd  part  of  |H(ju)|  will 
be  zero  If  |H(ju)|  Is  symmetrical  about  the  origin  ( 1 . e . ,  zero  clock- 
offset  and  aligned  antennas).  If  the  Input  process  to  the  filter 
H(Jw)  Is  symmetrical  about  the  origin  with  spectral  density  U(u) 

[constant  if  white  noise]  then  the  auto -spectrum,  SaaU)  and  cross-spectrum 
S  a(w)  of  the  output  quadrature  processes  o(t)  and  B(t)  ,  are  given  by 

ap 

[Frank,  1969] 

Saa(“)  "  S66<“)  *  Ev  |HU“)i2  U(-) 

?  (17) 
and  SaB(<-)  -  -j  00|H(j«)r  U(«) 

If  Input  complex  white  noise  is  correlated  and  is  of  asymmetric  spectrum 
(which  results  In  auto  spectral  density  Suu(w)  and  cross-spectral  density 
Suy(w)  of  the  input  quadrature  processes  u(t)  and  v(t)),  then  [Frank,  1969] 


-  Ev  H(joi)|2  Suu(o,)  +  j  O  |H(  jw)  I2  Svu  Coi) 
SaB(w)  =  -j  00  |H<»|2  Suu  (u>)  ♦  EjHfja,)!2  Sy(j  («) 


08) 


0. 


These  linear  transformations  on  auto-correlation  and  cross-correlation 
functions, not  on  the  signals  themselves, are  shown  symbolically  in  Fig.  3(a) 
and  (b).  Fig.  3(b)  results  from  3(a)  when  input  process  is  white 
uncorrelated  noise.  These  diagrams  make  it  clear  how  non-symmetrical  band¬ 
pass  characteristics  of  the  channal  Introduces  correlation  between  Inphase 

and  quadrature  components  of  the  band  pass  process.  Now  a  method  is 

o 

discussed  to  obtain  H(S)  from  |H(jw)|  and  the  eigen-values  of  H ( S ) 
are  then  related  to  the  tropo-link  parameters.  First,  a  general  complex 
model  will  be  derived,  and  this  model  is  then  applied  to  the  tropo-channel . 


Low-pass  equivalent  of  transformation  relating  input  and  output 
into  correlation  functions  when  input  Gaussian  noise  is  correlated 


Fig.  3(a) 
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R  (t) 
uu 
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Ev  |H<  jcu)  (2 

- > 

: —  i 

1 - > 

-j  00  J  H  (  ja>)  J 2 

- ^ 

R  (x) 


R  JO 


Low  pass  equivalent  of  transformation  when  Input  Gaussian  noise  Is 
uncorrelated 

Fig.  3(b) 


A  Complex  Model.  The  motivation  to  simulate  the  band-pass  channel  (filter) 
into  the  low-pass  equivalents  (or  complex  filter)  stems  from  the  fact  that 
the  inherent  low-pass  nature  and  speed  limitations  of  the  analog  computer 
make  its  operation  in  the  low-pass  mode  highly  desirable.  In  the 
following  discussion,  we  shall  show  how  a  band-pass  system  can  be  simulated 
equivalently  by  two  inter-connected  low-pass  systems. 

Consider  the  band-pass  filter  with  transfer  function  R(f)  and  its 
low-pass  equivalent  H(f).  The  input  signal  Is  given  by: 

x ^ ( t )  *  Re  [z'-j ( t)  exp(j2nfot)"| 

■(  1  /2 )  z  i  (t)  exp(j2JIf0t)  +(l/2)z]*(t)  exp  (-J2nf6t)  (19) 

where  z^(t)  =  u(t)  -  jv(t) 

Hence  X^f)  =(1/2)Z]  (f-fo)  +(l/2)Z1*(-f-fo)  (21) 

The  transfer  function  of  the  filter  can  be  expressed  In  terms  of 
H ( f )  in  a  similar  manner, 

R(f)  *(  l/2)H(f-fo)+(  l/2)H*(-f-fo)  (22) 


where 


H(f-fo)  =  (1+sgnf)  R(f) 


(23) 
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Now  for  the  output  signal, 

X2(f)  =  R(f)  X^f) 

={l/4)H(f-fo)  Z1  { f-fo)  +( 1/4)H*( f-fo)  Z^*(-f-fo) 
=(l/2)Z2(f-fo)  +(l/2)Z2*(-f-fo)  (24) 

From  eg.  (24),  It  Is  evident  that 
Z2(f)  =(l/2)H(f)  Z,(f) 

or  Z2(t )  =(  1  /2)h(t )  (XyZ,  (t )  (25) 

Output  *2(t)  In  complex  form  is  written  as 

Z2(t)  -  a(t)-j6(t)  (26) 


Now  we  decompose  h(t)  into  Its  real  and  imaginary  parts  so  that 


(l/2)h(t)  =  h^t)  +  j  h2(t). 


(27) 


where  h^t)  and  h2(t)  are  real  and  can  be  considered  as  inpulse  responses 
of  real  low-pass  filters.  The  corresponding  transfer  functions  are  given  by 


H,(f)  -  1/4  [h( f)  +  H*(-f)] 
and  H2(f)  »(l/4j)[H(f)  -  H*(-fj] 


(28) 


i.e.,  by  the  symmetrical  and  anti-symmetrical  parts,  respectively.  It 
follows  that  H-| (f)  and  Hg(f)  are  related  directly  to  the  parts  of 
the  band-pass  function  R(f)  which  are  symmetrical  and  anti-symmetrical, 
respectively  about  the  center  frequency  fo.  H^(f)  and  H2(f)  are 
related  to  even  and  odd  parts  of  |H(jw)|^  (eq;  16)  by: 


k 
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£v|H( j“) !'  *  |Hl(f) I2  +  |H2(f)j2 

and  0D|H(ju)|2  «  jh2(f)  H^f)  -  JH2*(f)  H^f)  (29) 

From  examples  (20),  (25),  (26)  and  (27),  one  writes: 
z2(t)  =  a(t)  je(t)  =  (u-jv)  ®  (h1  +  j  h2) 

=  (u®h)  +  v®h2)-j  [v  ®  h.  -u  X  h?] 
or  In  frequency  domain 

a(f)  =  U(f)  H1(f)  +  V(f)  H2(f) 

B(f)  =  V(f)  H1 ( f )  -  U(f)  H2(f) 

So  by  this  approach,  a  single  input,  single  output  band  pass  filtering 
operation  is  replaced  by  a  double  input,  double  output  low-pass  filtering 
operation  which  characterizes  separately  the  transfer  of  inphase  and  quadrature 
components.  The  equivalent  low-pass  system  is  represented  schematically  by 
the  lattice  arrangement  shown  in  Fig.  4. 


Fig.  4.  Equivalent  low-pass  filtering. 

•till' 

A  cascade  combination  of  band-pass  filters  can  be  handled  directly  In 
terms  of  low-pass  equivalents  as  shown  in  Fig.  (5).  [Frank,  1969] 


(30) 


(31) 


0.6i 


,(t) - 


Fig.  5  Low-pass  equivalents  of  cascade  combination  of  band-pass  filters 


We  should  note  that  if  a  band-pass  filter  exhibits  symmetry  about 
fQ,  then  Hg( f }  =  0  and  there  is  no  cross-coupling  between  the  In-phase 
and  quadrature  components.  This  Is  the  kind  of  transmission  required  to 
avoid  amplitude  to  phase  modulation  conversion  on  a  band  pass  signal. 

In  Appendix  1,  a  first  order  and  second  order  complex  system  is  derived 
which  when  driven  by  the  complex  white  noise  outputs  the  complex  narrow- 
band  process,  In  the  following  discussion,  the  eigen-values  of  the  complex 
system  are  estimated  from  the  first  three  spectral  moments  of  the  output 
Doppler  spectrum  (or  complex  spectrum).  First,  the  Individual  blocks 
[H-|-( j«) ,  HR(jw),  A(ju)]  will  be  simulated,  then  the  combined  system 
[H(jw)]»  simulating  the  tropo-llnk.wlll  be  synthesized. 

Simulation  of  the  Transmitter  Gain  Function  (equation  12).  Simulation  of 
the  transmitter  gain  function  Is  the  approximation  of  the  Gausslan- 
magnitude  approximation  given  by  equation  12.  Given  the  Ideal  Gaussian 
magnitude  characteristic : 


where  the  unit  bandwidth  has  been  selected  as  a  patter  of  convenience,  let 
find  an  all -pole  transfer  function 

H(s)  *  - ]— 

Gn(s,a) 

where  G  (s,7)  Is  a  real  Hurwitz  polynominal  with  a  coefficient  vector 

_  p 

a,  whose  spectrum  approximates  |Hy(jw)|  best  in  the  Taylor  sense  at  the 
origin, defining  the  error  function  as: 

E(s,a)  =  Gn(s,7)  Gn(-s,a)  -  exp[-s2] 

We  readily  see  that  the  optimum  G  (s,¥)  Gn(-s,7)  corresponds  to  the 
truncated  Maclaurin  expansion: 

A  6  2^ 

Gn{s,a)  Gn(-s,a)  =  l-s2  +  |r*|r+-*  +  (*l)n 

Finally,  Gn(s,7)  is  Isolated  from  the  product  Gn(s,a)  Gn(-s,a)  by 
factoring  and  retaining  only  LHP  factors.  Based  upon  the  above  analysis 
and  renormalizing  w.r.t.  unit  3-db  bandwidth,  Dlshal  [1959]  and 
Humphry  [1970]  obtained  the  pole  locations  up  to  the  order  n  =  10. 
Increasing  n,  obviously,  would  give  better  'matching'  characteristics. 

But  only  the  case  for  n  =  2  vtould  give  symmetrical  spectrum.  So  a 
simple  second-order  complex  model  Is  considered  and  It  Is  suggested  that 
further  work  is  necessary  to  obtain  better  'matching'  characteristics 
without  distorting  the  symmetry.  By  considering  the  pole  locations  at 
a  +  jb,  where  a  *  -  l*2841100aQ  and  b  *  *  j*531895a0  (oQ  Is  3-db 
bandwidth  given  by  equation 5),  one  obtains  the  second-order  complex 
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model  (see  Appendix  1 .equation  AT -7)  as  in  Fig.  6 


Uncorrelated 
complex  white 
noise 

(u1(t)+ju2(t)) 


l 

J 

(s+a)2  +  b2 

X2(t) 

'  l 

L . . 

^(t) 


(s+a)2  +  b2 
Fig.  6 


> 


Simulation  of  Angular  Dependence  (equation  14).  The  angular  dependence 
function  is  given  by  equation  14,  e.g.. 


where  K 

o 

obtains 


|ha(Ju)|2*  ■  °-  ■  2 

A  (f2  +  <2r 


is  given  by  equation  (7).  From  the  above  equation,  one 


H,(S) 


Z"v?/V^ 

(s  +  2nK0)2 


(32) 


The  simulation  of  H^( s )  In  quadrature  form  can  be  obtained  by  assuming 
a  *  c,  b  »  d  »  0  In  equation  Al-3  of  Appendix  1,  and  then  substituting 
In  equation  Al-7,  one  obtains: 
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H1 (s)  =  He(s)  and  H2(s)  =  0 


and 


u-,  ( t ) 

White  noise 
(u1(t)+ju2(t)) 
“2(t) 


HA(s) 

- 

- 

Ha(s) 

_ t) 


> 


BAU) 


Fig.  7 


Simulation  of  Receiver  Antenna  Gain  Function  (equation  13).  Two  methods 

are  discussed  to  simulate  the  receiver  antenna  gain  function  given  by 

equation  03).  It  Is  noted  that  the  transmitter  and  receiver  antennas  are 

Identical  except  that  the  receiver  antenna  is  not  pointed  at  the  great-circle 

but  is  offset  by  an  amount  a.  correspondina  to  a  Ooppler  offset  of  f  . 

a  a 

Method  1 .  In  this  method,  first  the  spectrum  |H^.(ju)|2  |HA(jw)|2 

(Gj(f)  A( f ) )  would  be  shifted  left  by  an  amount  fa.  The  shifted 

2 

spectrum  would,  then,  be  multiplied  by  |H^(ju)|  with  fa  »  0  (it  is 
equivalent  to  |Hj(jui)|  ).  Finally,  the  spectrum  would  be  shifted  right 
by  The  overall  simulation  Is  shown  in  Fig.  8.  The  spectral-shift 
unit  shown  in  Fig,  9  performs  a  single  side-band  modulation  directly  at 
baseband.  The  mathematics  of  the  signal  processing  Is  outlined  below. 

[Mans ke, 1968] 

Consider  the  case  In  which  a  narrow-band  signal,  expressed  In  the 
form  of  x(t)  ■  I ( t )  Cosu>ct+Q(t)  s1mwct  (33) 

is  passed  through  a  filter  having  characteristics  centered  at  a  frequency 
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«], which  Is  different  from  #»c  by  a  small  amount  (i^  -  uc  a  ffl).  The 
signal  can  then  be  expressed  in  terms  of  ,  rather  than  «c  in  the  form 
of: 


X{t)  =  apjOCosi^t  +  e^tjslnu^t 
From  equations  (33)  and  (34),  one  can  obtain  (see  Fig.  9) 


(34) 


Quadrature 

Inputs 


Left 

->}  Spectral 
Shift 

by  f 


(s+a)2  +  b^ 


1  (s+a)2  +  b2 


Right 
.^Spectral 
*  Shi  ft  | 

I 

by  f. 


_ ^  qr(  t) 

Quadrature 

outputs 


-BR(t) 


Fig.  8.  Simulation  of  receiving  antenna  gain  function  (method  1) 


a(t)  =  I(t)  Cos  ft  -  Q(t)  Sin  ft 

d  a 


and  6(t)  *  I(t)  Sin  f  t  +  Q(t)  Cos  ft 

a  o 


(35) 


<xKti) 


Fig.  9.  Spectral-shift  unit 


Method  2.  In  this  method,  the  receiver  antenna  gain  function  Is  simulated 
by  the  complex  model  as  described  In  Appendix  1.  If  In  equation  (Al-3), 


one  assumes  a  =  C  =  a ,  d  =  0,  b  ■  bo,  then  one  can  obtain  from  equation 


Al-7  that:  [See  Fig.  10] 
1 


Ms)  =  ,  9 

(s+a0)2+b2 


and  H.(s) 


2  (s+a0)[(s+a0)2+b02] 


(36) 


and  from  equation  (Al-11  and  12),  one  obtains 
b„ 


R(o) 


‘j  ~7 


and 


2. 


R(o)  -  -  (aS  -J-  ) 


(37) 


Since  the  first  moment  (gives  the  average  frequency  of  the  spectrum)  Is 
1  * 

equal  to  R(o)  and  the  second  central  moment  (related  to  f(o))  gives 

the  3-db  band-width,  one  derives  from  equation  37,  that 

b. 


-I 

fa  “  <-  3T  )  lr  or  b 


2n  w'  "o'*  4nfa 


(38) 


Similarly, 


0  (2n)2 


=  1 


(HT 


,4 


(39) 


From  equations  (38)  and  (39),  one  obtains 


*  (2n)z  [.*-  fa2_ 


(40) 


For  a*  >  o,  one  requires  that  |o0|:  >  |fa|  which  Is  Intuitively  obvious. 
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Fig.  10.  Simulation  of  receiver  antenna  gain  function  (method  2) 


Simulation  of  Clock-Offset  [COF '  operator.  Equation  1].  It  has  been 

mentioned  earlier  that  the  clock-offset  (  Af  )  merely  shifts  the  spectrum 

either  to  the  left  or  right  depending  upon  the  positive  or  negative 

clock-offset.  The  spectral-shift  unit  can  be  Implemented  similar  to 

that  shown  In  Fig.  9,  by  substituting  Af  Instead  of  f  .  Drift  In  the  clock 

« 

can  be  incorporated  by  varying  Af  accordingly. 


Over-all  Tropo-Channel  Transfer  Function  TpjffiKP  Scattering) 
The  overall  channel  transfer  function  In  quadrature  form  Is  the 
cascaded  simulation  of  various  terms  In  equation  1  and  Is  shown  In 
Fig.  11. 
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Fig.  11.  Over-all  tropo-channel  simulator  for  diffused  scattering 
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Using  Fig.  5,  one  can  obtain  the  cascade  form  of  the  transfer  functions 
of  transmitter  and  receiver  antenna  gain  functions  (Figs.  6  &  10) 
and  angular  dependence  function  (Fig.  7)  as  shown  In  Fig.  12,  where 


Hcr"  [(s+a)2+bZ][s+k3Z[(sM0)Z+b02] 

Hu"*  [(s+a)Z+bZ][s+k]Z(s+a0)[(s+ao)Z+boZ] 


(40) 


Fig.  12.  Cascade  form  of  over-all  tropo  channel  simulator  as 
shown  In  Fig.  11 . 

Models  for  Low  Frequency  Specular  Components: 

Duvoisin  (1969)  points  out  that  in  the  presence  of  multiple  scattering 
layers,  the  total  Doppler  power  spectrum  is  the  sum  of  the  spectral  contributions 
of  the  Individual  layers.  The  shape  of  the  portion  of  the  Doppler  (power)  • 
spectrum  of  the  received  signal  caused  by  a  given  horizontal  scattering  layer 
can  therefore  be  shown  to  be  controlled  by  the  large  numbers  of  scattering 
points  drifting  across  the  scattering  volume  and  Is  functionally  represented 
by  the  equation  2.  Considering  equations  (3  to  7),  Duvoisin  (1969)  concluded 


i 
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that  the  narrow  power-peaks  near  zero  Doppler  could  be  produced  by  any  one  of 
the  following  conditions: 

1.  Very  low  values  of  cross-path  wind,  V. 

2.  Very  narrow  beam  antennas  pointed  down  the  great  circle  azimuth. 

3.  Very  narrow  angular  dependence  function. 

Since  neither  conditions( 1 )  or  (2)  were  satisfied  during  most  of  the 
Barrow,  et  al.  (1968)  rake  sounding  experiments,  the  specular  peaks,  in  most 
cases,  were  probably  caused  by  condition  3.  The  angular  dependence  function 
used  in  Birkemeier,  et  a  1 . { 1 96 8)  to  simulate  the  narrow  specular  peaks 
observed  in  the  rake  data  was  of  the  form 


Vf>  = 


2n  A  t  1 

-P 

r,  k2i 

n  z 

f2  + 

7  j 

A2 

n  __ 

1  -P/2 


(41) 


where  kQ  is  given  by  eq.  7,  An  is  the  degree  of  ahisotropy, lz  is  the 
spatial  correlation  length  in  a  vertical  direction  and  p  is  ^  for  Tatar 
ski  scattering. 

For  An  =  1  (isotropic  case)  and  p  =  >  equation  41  reduces  to  the 

equation  6,  where 


Ao  = 


2n  t. 


-  11/3 


(42) 


l  J 

In  eq.  (41),  the  narrow-peaks  can  be  obtained  by  either  increasing  p 
or  increasing  An  .  According  to  Duvoisin  (1969),  maintaining  A  =  1  and 
increasing  p  could  produce  sufficiently  narrow  specular  peaks  only  when  p 
was  increased  by  almost  two  orders  of  magnitude.  However,  a  power  law  of 
p  =  11/3  and  a  degree  of  anisotropy  of  A  =  20  produced  peaks  having  the 
necessary  Doppler  bandwidth.  So  equation41  can  be  rewritten  as: 


Vf> 


*0 


2  k2„ 
f  +T 

A2 

n 


-11/6 


(43) 
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Similar  to  equation  8,  for  small  antenna  misalignment,  A(f)  can  be 
approximated  by: 


n(f)  =  A  [f2  +  k2  ]'2 
n  on  1  onJ 


(44) 


where  A  =  A  A 
on  on 


and  k  =  k  /A 
on  o  n 


(45) 


Now  similar  to  equation  32  and  Fig.  7,  one  can  obtain 
2n3/2/{A„  /k  ) 


HAn<s)  " 


(s+2nk  ) 

on 


on  on' 

r 


(46) 


and 


^(t)- 


1  HAn(s) 


o(t) 


“2(t) 


6(t) 


Fig.  13 


Now  similar  to  Fig.  11,  one  can  obtain  the  overall  tropo-channel  simulator 
for  low  frequency  specular  components  i.e., 


U,(t)_  _ 

Uncorrelated 
white  noise 

u2(t)  - 


Fig.  8 

, 

Fig.  6 

Fig.  13 

or 

clock-offset 

Transmi tter 

Angul ar 

Fig.  10 

/ 

operator 

antenna 

dependence 

°eceiving 

Pi  n  Q 

gai  n 

functi on 

antenna  gain 

— ■ • 

function 

_ _ J. 

function 

> 

-> 


-V 


a 


f 


Fig.  14.  Overall  tropo-channel  simulator  for  low-frequency 
specular  components. 


Overall  tropo-channel  transfer  function  can  be  obtained  by  combining 
Fig.  (11)  and  (14)  and  is  given  in  Fig.  15. 


11  ">  I.I,SS--|'|_.J|  LinUBIrJBSS 
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Model  for  the  Average  Drift  of  the  Scatters: 

In  the  simulation  of  the  angular  dependence  function,  symmetrical  form 
of  the  function  is  assumed.  Since  we  would  expect  on  the  average  as  many 
up-Doppler  signals  as  down-Doppler  signals  and  therefore,  the  average  total 
doppler-shi ft  for  aligned  antennas  would  be  zero.  If  there  are  more  scatters 
on  one  side  of  the  great  circle,  it  results  in  the  average  frequency  (f^)  to 
be  non-zero  and  the  simulation  of  the  angular  dependence  function  (as  shown 
in  Fig.  7  and  13)  would  be  modified.  Using  the  results  of  the  Appendix  1 
(Al-15) ,  one  obtains : 


a^(  t ) 


BA(t) 


Fig.  16 


where  H^(s)  = 


z”3/2A7 


on 


1 

(s+2nkQn)2+AA2 


Lw/ 

i  / 


A„  /k 
on  on 


6A 

(s+2nkon)2+AAS 


(46) 


H2(s)  = 


Now  similar  to  equation  (38),  one  can  obtain: 

aA  =  2nfd  (47) 

If  the  horizontal  drift  of  the  scattering  is  occuring  in  a  deterministic 
or  random  fashion,  tien  aA  (the  imaginary  part  of  the  complex  pole 
-ic  =  2nkQn  +  jaA)  can  be  varied  accordingly. 

Examples : 

In  Fig.  12  a  tropo-scatter  channel  simulator  for  diffused  layer  is 
shown.  The  poles  of  this  channel  transfer  function  are  related  to  the  tropo- 
link  parameters.  The  variations  in  the  pole  locations  would  correspond  to 
various  link  conditions,  resulting  in  different  spectral  moments  as  described 
in  earlier  reports.  In  the  following  ,  for  illustrative  purposes,  the  pole 
locations  [-k^.-k.,]  of  a  2nd-order  complex  system  are  varied  to  simulate 
various  link  conditions  and  their  effect  on  the  spectral  properties  of  the 
output  quadrature  processes  are  tabulated  in  Table  1  .  The  Doppler  spectrum 
(eq.  Al-10)  and  auto  and  cross-correlation  functions  (real  and  imaginary  parts 
of  the  equation  Al-9)  are  plotted  in  Fig.  17  and  18  for  all  the  cases  described 
in  Table  1.  Some  of  the  conclusions  from  this  analysis  are  summarized  below. 

1.  For  symmetrical  spectrum  with  no  clock-offset,  the  cross -correlation 
function  is  zero  for  all  t  (case  1,2, 7 ,8). 

2.  For  symmetrical  spectrum,  the  clock-offset  (average  freq.)  is  obtained  by 
averaging  the  imaginary  part  of  the  complex  poles  (case  3, 3c, 6, 6c). 

3.  Slope  of  the  cross-correlation  function  at  zero  by  [Rag(°)]  is  obtained 
from  the  equation  Al-11,  and  is  equivalent  to  the  average  frequency  of 
the  spectrum,  so  the  sign  of  Rag(o)  can  be  computed  from  the  average 
frequency  of  the  process  (1st  moment). 


For  asymmetrical  spectrum,  the  slope,  ft  „(o),  could  be  zero.  This 

ap 

results  due  to  the  fact  that  the  average  frequency  of  the  asymmetrical 
spectrum  due  to  the  'asymmetry'  alone  is  equal  and  opposite  in  sign 
of  the  clock-offset.  So  the  net  average  frequency  of  the  spectrum  is  the 
sum  of  these  two  effects  and  is  therefore  zero  (case  4, 4c, 9). 

The  auto  correlation  function  of  the  quadrature  processes  after  hard- 
clipping  decays  faster  than  the  case  of  the  unclipped  signals. 
Hard-clipping  of  the  processes  does  not  change  the  zero-crossings  in  auto 
and  cross-correlation  function. 


D.  75 


Spectrum 

symmetrical 


Case 

no. 

Pole  locations 

(-VV 

& 

asymmetrical 

Average 

frequency 

Cross-correlation 
P>  „(t)  for  all  i 

dp 

Sense  of 

1 

-1+j  ,-1-j 

symmetrical 

zero 

zero 

zero 

2 

-l+2j,-l-2j 

symmetrical 

zero 

zero 

zero 

3 

-1  ,-1-j 

symmetrical 

(-ve) 

non-zero 

(-)ve 

3c 

-1,-1+j 

symmetrical 

( +)  ve 

non-zero 

(+)ve 

4 

-2+2 j  ,-1-j 

asymmetrical 

zero 

non-zero 

zero 

4c 

-2-2 j  ,-1+j 

asymmetrical 

zero 

non-zero 

zero 

5 

-2  ,-1-j 

asymmetrical 

(-)ve 

non-zero 

(-)ve 

5c 

-2, -1+j 

asymmetrical 

(  +  )ve 

non-zero 

(  +  )ve 

6 

-2,-2-j 

symmetrical 

(-)ve 

non-zero 

(-)ve 

6c 

-2 ,-2+j 

symmetrical 

(  +  )ve 

non-zero 

(+)ve 

7 

-2+j  ,-2-j 

symmetrical 

zero 

zero 

zero 

8 

-2+2 j ,-2-2j 

symmetrical 

zero 

zero 

zero 

9 

- 3+3 j ,-2-2j 

asymmetrical 

zero 

non-zero 

zero 

10 

~  3+2 j ,-2-j 

asymmetrical 

(  +  )ve 

non-zero 

(  +  )ve 

Table  1.  Summary  of  the  Spectral  Analysis  of  Output  of  the  Complex 
System  for  Different  Pole  Locations. 
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Appendix  1 


In  this  Appendix,  a  new  concept  of  complex  modeling  is  introduced. 

This  complex  modeling  (for  complex  processes  or  orthogonal  processes)  is 
shown  to  be  a  special  form  of  modeling  bivariate  processes.  The  complex 

representation  of  narrow-band  processes  has  already  been  described  in 

earlier  reports.  Here  the  state  variable  description  of  these  processes 

is  considered  and  its  application  to  modeling  the  scatter  portion  of  the 

channel  is  described. 

In  this  complex  state  variable  description  [Baggeroer,  1970],  a  complex 
random  process  is  generated  by  exciting  a  linear  system  having  a  complex 
state  variable  description  with  a  complex  white  noise.  The  only  assumptions 
made  are  on  u(t)  and  x(TQ)  where  u(t)  is  the  complex  white  noise  and  x(TQ) 
is  the  complex  initial  conditions  of  the  linear  system  which  could  be 
assumed  to  have  the  following  state  variable  description. 


x(t)  =  F(t)x(t)  +  G (t)u(t)  [Linear  State  Equation] 

(Al-1) 

z(t)  =  c(t)x(t)  [Linear  Observation  Equation] 
and  z(t)  =  a(t)  +  jg(t) 

where  all  the  coefficient  matrices  may  be  complex. 

With  the  assumptions  on  u(t)  and  x(T0),  one  can  obtain  the  covariance 
function  kj(T)  of  the  outPut  [Baggeroer,  1970].  From  the 

complex  spectrum  of  z(t),  s2(w),  can  be  obtained.  It  can  be  shown  that  the 
complex  output  z(t)  *  a(t)  +  jfi(t)  satisfies  the  conditions  of  narrow-band 
processes  (orthogonal  processes). 


In  this  way,  a  complex  process  can  be  generated.  It  will  be  shown 
below  that  the  different  coefficient  matrices  could  result  In  different 
complex  models  which  will  mimic  the  different  channel  conditions  resulting 
due  to  clock-offset,  misaligned  antennas,  varying  fade  rate  and  spectral 
shape  of  the  channel.  In  the  description  below,  an  example  (2nd  order 
system)  Illustrates  how  these  conditions  of  the  channel  can  be 
simulated.  This  example  is  described  In  [Baggeroer,  1970],  but  its  appli¬ 
cation  in  channel  modeling  is  pointed  out  here.  Finally,  lst-order  system 
is  considered  which  would  be  useful  in  simulating  angular  dependence 

function. 


Derivation  of  the  Channel  Transfer  Function 

Let  us  consider  a  second-order  system.  In  the  steady-state,  it 
corresponds  to  a  system  with  two  poles.  The  pole  locations  -k1  and  -k2 
need  not  be  complex  conjugates  of  one  another.  For  this  system,  the  state 
and  observation  equations  are: 


d_ 

dt 


and  z(t)  3  [1:0] 


X^t) 

0  1 

X2(t) 

“^1^2  ~ Ik] 2 ) 

r 

Xj(t) 

0 

+ 

x2(t) 

1 

u(t) 


x^t) 


x,(t) 


By  simple  manipulations,  it  can  be  shown  that  the  above  equation 
corresponds  to 
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Let  z(t)  =  alt)  +  j8 ( t ) 
=  a+jb 
k2  =  c+jd 
A  =  ac-bd 
8  =  bc+ad 
E  -  a+c 
F  =  b+d 


( AT  - 3) 


The  above  equation  ( A 1-2)  can  be  rewritten  as: 

alt)  +  jS(t)  ■  -(E+jF) (a(t)  +  jB(t)]  - 

[A+jB]-[a(t)  +  j@( t) ]  +  U](t)+jU2lt) 


(AT -4) 


Separating  real  and  imaginary  parts  in  Eq.  (Al-4),  one  can  obtain 


a(t)  +  Ea( t)  -  F6(t)  ♦  Aa(t)  -  Belt)  =  11,(0 
Bit)  +  EB ( t )  +  Fi( t)  +  Balt)  ♦  ABlt)  =  U2(t) 


(A1  -5) 


Equation  (Al-4)  can  be  written  in  transfer  function  form.  Let  us  assume 
that  Uj(t),  U2(t),  a(t).  Bit)  are  general  input-output  of  a  bivariate  system, 
such  that  their  Laplace  Transforms  can  be  taken.  Then  by  simple  manipulations, 
one  can  show  that 


als) 


(s^ES+A) 


(s2+ES+A)2*{FS*B) 


7  Ul(s) 


F  ^ - 2  U2(s) 

(s2+ES+A)Z+(FS+B)z 


and 


e(s) 


s2+SS+A 

(s2+ES*A)2+(FS+B)2 


U2(s)  ' 


-5-^ - y  u,ls) 

(sZ+ES+A)Z+(FS+B)<: 


(A  1-6) 


0.81 


s2+£S+A 


If  H,(s)  =  -5-  .  and  H.(s)  =  -y  - , 

<s  +ES+AJ  ^+(FS+B) ^  2  (s2+ES+A)2+(FS+B)2 


then,  equation  (Al-6)  can  be  written  as: 


(A1-7) 


c(s)  =  «1 (s)U, (s)  +  H2(s)U2(s) 
and  e(s)  *  H1(s)U2(s)  -  H2(s)U](sJ 

Equation  (A  1- 8)  can  be  written  in  block-diagram  form  as: 


(A  1-8) 


In  the  above  block  diagram,  it  can  be  observed  that  the  forward  transfer 
function  Hj(s)  in  both  the  paths  are  the  same  and  the  transfer  function  in 
the  cross-paths  [H2(s)J  are  equal  but  of  opposite  sign.  The  author  points 
out  here  that  this  Is  the  condition  for  modeling  complex  processes  (orthogonal 
processes)  in  quadrature  form. 

If,  in  equations  (Al-6)  and  (A1-8),  B*0  and  F=0,  then  H2(s)  is  zero.  For 
this  condition  there  is  no  cross-coupling  between  two  quadrature  components 
o(t)  and  e(t),  and  these  two  processes  are  uncorrelated.  As  discussed 
earlier,  zero  correlation  Implies  symmetrical  spectrum  with  no  clock- 
offset.  If  B=0  and  F=0,  then  a=c  and  b=-d  e.g.,  both  the  poles  are  mirror 
images  of  one  another  about  o-axis  in  s(=o+ju)  plane. 
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For  non-zero,  B  and  F  constants,  H2(s)  is  not  zero  and  a(t)  and  e(t) 
processes  are  correlated.  Since  B  and  F  are  functions  of  the  poles  of  the 
system  (. .kg)  >  that  implies  that  the  pole  position  determines  the  shape  of 
the  power  spectral  density  of  the  complex  process.  Equation  (Al-E)  can  easily 
be  simulated  on  an  analog  computer.  By  varying  pole  locations  (k^.k^), 
one  can  obtain  different  spectrums.  Also  note  that  the  differential 
equations  involving  a(t)  and  g(t)  in  Eq.  (Al-5)  are  coupled  equations.  So 
'coupled'  systems  (e.g.  non-ideal  transformer)  can  be  analyzed  by  this 
complex  or  phasor  analysis. 

The  covariance  function  K~(t,t)  of  the  output  is  given  by 
[Baggeroer,  1970]:  (  Pis  the  input  noise  power). 

(  k*  ek1T-k*  ek2T) 

P 

"(K2*kl)#  +  j.  J^lM  (eklT  -  ek2T 
Retki+i^] 

(i2  rt') 

P 

(kVkl}  Im[k1k23  [e"klT  -  e'k2T 

Re[k|k2] 

From  e.g.  (A-9),  one  can  compute  s^(w)  as: 

k,+kT  Re[k,] 

Let  A,  =  -4— ~  ’  - r—r~ 

(kg-^)  Refk^kg] 

k,+k»  Re[k, ] 

A  *  -  ■  f  *  - J-;- 

(k2-k.j)  ReCk-j+kg  3 


D. 


then 


S;(w)  =  2P 


Re[A-|]Re[k^]+Im[A^  ](u+Im[k^l) 


Re  ]+[u)+Im(k^  )]^ 


ReCAglReCkgl+Im^Kw+Imfi^]) 
Re  C^+L  +  Im^)]2 


(Al-10) 


From  e.g.  (Al-9)  one  can  obtain: 


r(0)  =  -  j  [Jig*!]  --  j| 


(A1-11) 


and 


*r*(0)  =  -  -a.(.?2tj2.).!£(a2>b2) 
'  '  a+c 


(‘  *  ¥) 


(Al-12) 


First-order  Complex  System: 

The  equations  that  describe  the  Ist-order  system  with  one  pole 
(-k)  are: 


and 


dx(t)  =  -  k  x(t)  +  u(t) 
dt 

z( t)  =  x(t) 


(state  equation) 


(observation  equation) 


(Al-13) 


Following  the  steps  similar  to  2nd-order  system,  one  can  obtain 
by  letting 


S  =  a+jb 

and  z(t)  ■  a(t)  +  Je(t) 


(Al-14) 
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that 


H,(s) 


s»a 

(s+a)Z+b^ 


and 


H2(s) 


(s+a)2+bZ 


(Al-15) 


Similar  to  equations  (Al-9)  and  (Al-10)  one  obtains  [Baggeroer,  1970] 
k-(T)  =  p  e*Retk]|x|  e-jIm[k]T 


(Al  -16) 


and 


2Re[k]P 

(u-Im[k])2+Re[k]2 


(Al-17) 


The  auto  and  cross  spectra  of  a(t)  and  e ( t )  can  easily  be  found  In 
terms  of  even  and  odd  parts  of  s^(u>).  Note  that  the  net  effect  of 
allowing  a  complex  gain  in  e.g.  (Al-13)  is  that  the  spectrum  has  a  frequency 
shift  equal  to  the  imaginary  part  of  the  gain. 
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Fig.  17a 
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Fig.  17b 
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Fig.  17c 
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Fig.  17d 
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Fig.  18d 
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